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Abstract 

Oh' The presence of Kaluza-Klein particles in the universe is a potential manifestation of string theory 



cosmology. In general, they can be present in the high temperature bath of the early universe. In 



CN ' particular examples, string theory inflation often ends with brane-antibrane annihilation followed 

OO , by the energy cascading through massive closed string loops to KK modes which then decay into 

in 

I lighter standard model particles. However, massive KK modes in the early universe may become 

■ dangerous cosmological relics if the inner manifold contains warped throat(s) with approximate 

|: _.,,...._...„„,„.™c„...,,.,..„,-^ 

massive glueballs of various spins become the dangerous cosmological relics. The decay of these 
^ ' angular KK modes/glueballs, located around the tip of the throat, is caused by isometry breaking 

which results from gluing the throat to the compact CY manifold. We address the problem of these 
angular KK particles/glueballs, studying their interactions and decay channels, from the theory 
side, and the resulting cosmological constraints on the warped compactification parameters, from 
the phenomenology side. The abundance and decay time of the long-lived non-relativistic angular 
KK modes depend strongly on the parameters of the warped geometry, so that observational 
constraints rule out a significant fraction of the parameter space. In particular, the coupling of the 
angular KK particles can be weaker than gravitational. 
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I. INTRODUCTION 



Often when we embed a new layer of elementary particle theory in the theory of the early 
universe, we find that new particles become dangerous cosmological relics. The study of this 
problem is one of the few tools available for theorists to bridge these two theories. For exam- 
ple, the GUT theory incorporated in the Hot Big Bang model predicted the overproduction 
of monopoles; SUGRA theory predicted the overproduction of gravitinos and moduli fields. 
It is interesting to explore if the embedding of string theory in cosmology results in new 
dangerous rehcs. 

Kaluza- Klein (KK) particles in the universe is a possible signature of string theory physics. 
Relic KK particles may be good or bad for cosmology, depending on their properties. If their 
abundance multiplied by their mass is tuned to be right, and their life time is significantly 
longer than the age of the universe, the KK particles can be dark matter candidates. Oth- 
erwise, and more typically, they are dangerous relics, which may overdose the universe, 
or violate Big Bang Nucleosynthesis or other astrophysical constraints, depending on their 
abundance and lifetime. Therefore, the theory of relic KK particles is a powerful tool to 
study string theory cosmology and to confront it with observations. 

The AdS/CFT duality allows us to translate from the language of the gravitational KK 
sector to gauge theories. Therefore, independently of string theory, the story of dangerous 
cosmological relics repeats itself in the context of gauge theories with extra isometrics, which 
contain massive glueballs of different spins. 

We will consider KK particles in a string theory and string theory cosmology setting, 
which drew significant interest in the last few years, and which can be advanced quantita- 
tively. Indeed, important progress has been made in the studying type IIB string theory 
compactifications with fluxes where the moduli are stabilized l|, These developments 
in string theory have raised interests in the embedding of early universe inflation and low 
energy physics in a fundamental and phenomenologically viable framework. 

Constructions with fluxes and conifolds induce warped regions of spacetime 
- throats - which in turn may be used to generate mass hierarchies. For instance, the 
hierarchy between the string scale and the electroweak scale may be generated a la Randall- 
Sundrum [5] when the Standard Model degrees of freedom are localized at the bottom of a 
sufficiently warped throat. Furthermore, warping is an essential ingredient in the attempts 



2 



to realize brane inflation in the slow- roll regime |6| or the DBI regime [7| , or in the fast-roll 



regime 



see also e.g. [9|, [id, for recent works. 
Warped brane/anti-brane inflation involves a pair of brane/anti-brane D3 and D3, which 
are string theory objects: although inflation in this model may be described with an effective 
four dimensional field theory, its end point, reheating, is associated with the D3 — D3 annihi- 
lation, which is not captured by field theory. Reheating from brane/anti-brane annihilation 
is rather different from (p)reheating in quantum field theory 1^, and was investigated in 



several papers [13| 



14 



15 



16 



17l |. which identified the channels of inflaton energy decay. 



cascading from tachyon annihilation to the lighter standard model particles through massive 
closed string loops, KK modes, and brane displacement moduli. Indeed, if the internal space 
has warped regions (throats), the mass of the KK modes associated with them is much lower 
than the fundamental scale, being redshifted by the warp factor at the bottom of the throat. 
If the KK modes are sufficiently light then it is natural to expect them to be produced during 
reheating. For instance, in models of brane/anti-brane inflation, KK modes are copiously 
produced by the decay of the closed string loops from the brane/anti-brane annihilation at 
the end of inflation, and their energy must be efficiently converted into Standard Model 
(SM) degrees of freedom. 

Notably, reheating is then bottlenecked through KK modes which are localized around 
the tip of the throat. Warped brane inflation and reheating after brane annihilation were 
considered for throat models with high isometries, often described by the Klebanov-Strassler 
(KS) solution j^. The KS geometry has S'0(4) x Z2 isometries around the tip of the throat, 
and it is close to AdS^ x T^'^ away from the tip. Other known examples of throat solutions 
also have similar isometries 18|. To make the inner space compact, the KS throat is glued 



to the unwarped Calabi-Yau (CY) part of the manifold. For cosmologists, a significant 
advantage of this model is that the metric of the throat is known analytically, so that KK 
modes can be studied in some detail. However, as it was found in |14|, the advantage in 
mathematics is in this case a disadvantage in physics. Indeed, the KK modes in internal 
spaces with isometries have conserved angular momentum quantum numbers associated with 
these isometries. The final decay products of the KK particles, the SM particles, do not 
have such quantum numbers. Therefore KK modes with (internal) angular momenta cannot 
decay into SM particles. Moreover, KK particles with opposite angular momenta cannot 
annihilate completely into SM particles due to the expansion of the universe. Therefore, we 
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encounter the severe problem that the freezed out fraction of the long-hved angular Kaluza- 
Klein modes can overdose the universe [l^ ^! The problem of angular KK modes from 
an inner geometry with throats is not specific to the warped brane/anti-brane inflationary 
scenario. Independently on the particular realization of inflation, light KK modes may be 
produced during reheating if they are sufficiently coupled to inflation. In any context where 
the reheat temperature of the very early universe is large enough, massive KK particles can 
be produced in the thermal bath. Then, KK modes with angular momenta will again freeze 
out, and can be dangerous for cosmology. 

In a sense, the problem of dangerous angular KK modes in string theory cosmology is 
a re-appearance of the problem of KK modes from dimensional reduction in supegravity 
in higher dimensions with isometrics [20i]. However, the string theory setting brings new 
features in cosmology with KK modes. The unwarped CY part of the compact manifold 
does not preserve the isometrics of the throat. Therefore, when we embed the KS throat into 
the compact manifold, the isometrics of the throat inevitably get broken. These deviations 
from isometrics are significant at the base of the throat, close to the unwarped CY, but they 
decrease exponentially towards the tip of the throat. On the other hand, the KK modes are 
localized at the tip of the throat and their wave functions decrease exponentially towards 
the base. The net effect will be some violation of the conservation of the angular momentum 
of the KK modes. In terms of the four-dimensional effective lagrangian, this manifests itself 
by additional terms describing the decay of angular KK modes into SM particles with small 
coupling constants, controlled by the isometry breaking. This makes angular KK modes in 
the universe not stable but short or long lived, depending on the parameters of the model. 

In this paper, we investigate in some detail the decay of Kaluza-Klein modes with angular 
momenta associated with the compact manifold in warped flux compactifications, namely, 
in a KS throat with broken isometrics. In Section [TTl we begin with the background throat 
geometry; we consider first the KS solution with intact isometrics, and we then introduce 
the isometry breaking. We use the static perturbations of the KS metric considered in 21| . 



where they were connected to symmetry breaking operators in the context of the AdS/CFT 



^ It turns out that the left-over abundance of angular KK modes is very sensitive to the parameters of the 
throat geometry. Alternatively, after tuning the parameters, these modes have been proposed as possible 
dark matter candidates 
other models). 



14| . 15| (see also 19| for earlier proposals of KK dark matter, in the context of 
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correspondence. 

In Section IIIII we consider the wave functions of the KK modes in the throat geometry. 
For definiteness we will discuss the gravi-tensor KK modes (the spin-two KK modes, or 
spin-two glueballs). The knowledge of the wave functions allows us to integrate the internal 
space out, and to study the effective four-dimensional physics of the model. From the 
spectroscopy of the "KK atoms" (i.e. their masses as a function of the quantum numbers 
associated with the internal dimensions) we can constrain the interactions of the KK modes 
with each others and with SM particles; this results in "selection rules" which control the 
decay of these modes. We focus on the AdS part of the geometry as an approximation of 
the actual KS throat. We employ standard (quantum mechanics) perturbation theory to 
compute the effect of the isometry breaking perturbations on the wave functions of the KK 
modes. 

The interactions of the KK modes are studied in Section IIV[ We consider the model 
where the SM fields are localized on a 3-brane around the tip of the same throat where 
the excited KK modes are located. We study the possible decay channels of the long-lived 
modes and find that the leading one is into SM particles. Note that this channel differs 
from the deca y o f scalar KK modes into massless gravi-vector particles which was recently 
considered in 27l |. 



The second part of our paper is devoted to the consequences that the KK modes can have 
for cosmology and astrophysics. We first, in Section |Vl estimate the abundance of angular 
KK modes which freeze out in the hot expanding universe, and also estimate their life-times. 
We will mostly work with the conservative assumption of initial thermal equilibrium of KK 
modes. 

It turns out that the theoretical predictions for the angular KK modes abundance and life- 
time are very strongly dependent on the parameters of the modes. Therefore, in Section |VT] 
we put together observational constraints on the parameters of KK modes and consequently, 
on the underlying throat model, using limits from BBN 29|] and the astrophysical 7-ray 
background jSlJ, in case of decaying modes, and limit on the dark matter energy density, in 
case of lifetimes much greater than the present age of the universe. 

We restrict our calculations of the angular KK modes to single-throat settings. There 
are several reasons for that. The case of two throats, with one of them associated with 
(brane-antibrane) inflation (of energy scale 10^^ GeV or so) and the other one associated 
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with the standard model sector (of TeV energy scale), is often considered in the literature. 
In the cosmological context, however, the story of KK modes in the two throat models with 
such significant energy offset has the following conceptual complication, which usually is 
not addressed in the literature. Indeed, first, the Hubble parameter scale H ~ 10^'^GeV of 
infiation exceeds the TeV mass of the SM throat KK modes. Similarly, the temperature T 
after brane-antibrane reheating exceeds that low mass scale. For those situations the low 
energy SM throat will be screened by a horizon ^. In other words, the KK mass spectrum of 
the low energy throat has a gap below defined by H or T. Therefore, the simple picture of 
KK modes tunneling from the infiationary throat to the SM throat throat is relevant only 
after the temperature of the universe drops below the TeV scale, i.e. when the timing of 
tunneling is longer than 10~^^ sec. Another reason is that tunneling of the angular KK modes 
between the throats through the bulk of the inner manifold is exponentially suppressed. 

In the Discussion Section we summarize our calculations and discuss generalizations of 
our results to other SM models and multiple-throat cases. 



II. BACKGROUND WARPED GEOMETRY 



In this section, we discuss the background geometry with the warped throat as a part of 
the inner space manifold. We will use this construction in the rest of the paper. 

We consider a type IIB compactification with fiuxes turned along the internal dimensions. 
A common feature of fiux compactification is warping. The 10-D metric can be written 



or alternatively 



ds"" = e^^^y^ g^, dx'^dx" + gM dy'^dy" . (2) 



The greek indices /x, z/. A, ... run over the 4-dimensional coordinates of the external space- 
time, while the latin indices a, 6, c, ... run over the 6-dimensional coordinates y"' of the internal 
compact space (whose metric is cjab)- The outer space of the 10-dimensional background is 
approximated with the 4-dimensional Minkowski metric, gf^i, = rj^^j . 



^ The emergence of the horizon in the KS throat was considered, e.g. in 3C|. 
^ This approximation is poor when the Hubble parameter of the expanding universe exceeds the typical 
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The 10-dimensional geometry involves one or several "throats", each of them on one 
side ends with a regular "tip", and on the other side is smoothly glued to the bulk of the 
CY manifold ^. The known solutions for the throat, e.g. the KS solution, involve specific 
isometries and formally have a non-compact radial dimension. In order to get a sensible 
four-dimensional effective theory, the throat has to be glued to a compact CY which does 
not preserve these isometries. This results in deviations from the isometric KS solution, 
which are increasing as we approach the bulk of the CY. On the other hand, these isometry 
breaking perturbations are expected to be exponentially suppressed at the bottom of the 
throat. 

In the following sub-section, we discuss in more details the geometry of the throat. In 
sub-section III B[ we consider the breaking of the isometries due to the compact CY. In 
sub-section III C\ we discuss the essential parameters of the model. 

A. Isometric KS throat 

The known example of analytical throat solution ending with a regular tip is given by 
the Klebanov-Strassler solution [sl], where the warped throat contains a deformed conifold. 
The deformed conifold ^] is a 6-dimensional space defined in = M.^ by the quadric 

(3) 

1=1 

where Zi are four complex coordinates. Far from the origin Zi = 0, the deformation e may 
be neglected and the space approaches a cone. Its base, given by the intersection of ([3]) with 
the sphere in C^, 

4 

\zi\'^ = constant , (4) 

1=1 

is the 5-dimensional manifold T^'^, a fibration of U{1) over S"^ x S"^, which has the topology 
of ^2 X S^. 

mass scale at the tip of the throat (hke the masses of KK modes). In this case the tip of the throat 
would be screened from the rest of the compact manifold by an horizon and the KK mass spectrum would 
acquire a mass gap >, H. In this paper H <^ m ^ T . 
^ The embedding geometry of the inner manifold reminds us of "an udder" with multiple tips sticking out. 
Baroque "udder" constructions are not so baroque for relativistic astrophysics: recall that the embedding 
diagram of the multiple supermassive black holes in the galaxies around us in an expanding universe also 
looks like an udder! 
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The warping is achieved by turning on M units of RR flux through a 3-cycle of the 
deformed conifold and —K units of NS-NS flux through the dual 3-cycle. The full KS 
solution sl involves rather tedious formula for the metric. Not all of its information will be 
needed for our discussion, therefore we will work with simpler approximations. Away from 
the bottom of the throat, the geometry is approximately given by the Klebanov-Tseytlin 



solution {4] 



where dsf is the line element of T^'^. The warping may be written 



Rl + Rt In 



Hir) = -^^^ (6) 

with 

Rl = ^a''gsMK and Rt = — a'"" M\ (7) 
4 8 

where y/a' is the string length, gg is the string coupling, and the amounts of fluxes M and 
K are chosen such that -R+ > R^ (see Eq. (11 01) below). This approximate solution is valid 
only in the range R^ e~^+^^-^^ < r < R^. At smaller r, the singular metric based on 
the conifold has to be replaced by the regular one based on the deformed conifold jsj]. At 
r ~ _R+, the geometry is glued to the compact CY. 

Locally, the metric ([5]) corresponds to the direct product of AdS^ and T^'^ with the same 
radius of curvature R, but with R'^ = R^ + Rt In j slowly (logarithmically) varying 
from R = i?_ to R = R+. Taking i? to be a constant between i?_ and the metric 
reduces to the simple form 

ds^ = r/^^ dx^'dx" + dy^ + R^ fijiVL^) dO'dO^ , (8) 

where the coordinate y = R \n{R/r) measures the proper distance in the radial direction, and 
goes from y = to y = yt- fij is the metric on a 5-dimensional manifold X^, approximated 
by T^'^ in the present case, and i,j,k,... run over the 5 "angular" coordinates 9^ (denoted 
collectively as ^5). 

At the bottom of the KS throat (corresponding to y > yt in (^), the T^'^ evolves into 
a round of finite radius, while the two other angular directions, corresponding to the S'^ 
fibered over the S"^, shrink to zero size. The geometry may be described in terms of a new 
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radial coordinate p, with r oc e^/^ at large p. The asymptotic behavior of the metric in the 
vicinity of the tip at p = reads as 

where the warp factor e~^^'/'^ and the radius of are approximately constant, and 
eij{6^) denotes the metric of the S'^ fibered over the S^. The total variation of the warp 
factor from the top to the bottom of the throat is approximately given by ll] 

which can generate naturally large hierarchies for suitable choices of fluxes. 

The isometry group of T^'^ is SU{2) x SU{2) x f/(l), as may be seen from its embedding 
dl]) and (El) with e = 0, where SU{2) x SU{2) ~ 50(4) rotates the Zi's and U{1) multiply 
them by a phase. This is broken into SO{A) x Z2 by the deformation e of the conifold, where 
SO {A) corresponds to rotations on the 5''^. These are the isometrics of the KS throat. 



B. Breaking Isometries by Gluing to CY 

Kaluza-Klein modes in the 4 + 6 dimensional KS geometry with isometries in the angu- 
lar dimensions carry angular momenta quantum numbers associated with them. Standard 
Model particles in the 4-dimensional outer spacetime do not have those quantum numbers. 
Therefore in an expanding universe excited angular KK modes will never be converted com- 
pletely into the SM particles. In the model with isometric throat angular KK particles are 
the dangerous cosmological relics. 

However, when the throat is glued to the compact CY at y = 0, its geometry deviates 
more and more from ([8]) as ?/ decreases towards zero. In this sub-section we consider the 
background geometry of the throat which is smoothly glued to the big compact manifiold. 
We will model this geometry as static perturbations around the isometric geometry ([8]). We 
write ^ 

ds'' = e-2J'/^ [1 + e{y) w{n,)] r,,, dx^dx'' + dy^ + [f,,{Q,) + e{y) Sfi,{Q,)] dOW , (11) 

^ In Eq. (fTTjl . the perturbations e{y)w{fl5) and e{y) S fij (fl^) should be understood as sums over pertur- 
bations with different radial profiles, see Eq. ([T2)) and the discussion below. We will sometimes keep the 
"symbolic form" e(jj)w{fl5) and e{y) S fij (fl^) to simplify notations. 
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where SfijlQ^) and w{fl5) are the isometry breaking perturbations at a given y of the metric 
on X5 and of the warp factor respectively. For generic perturbations, only the 4- dimensional 
Lorentz invariance remains in fllip . The perturbations are expected to be of order one at 
the base of the throat, where it is glued to the bulk of the CY, and to decrease exponentially 
away from it, e{y) ~ q-^'v/^. In general, other supergravity fields and other components of 
the 10-dimensional metric may be perturbed as well, however, the form (fTT]) will capture 
the effects we are interested in. 

In principle, for a given throat model, the form of the isometry breaking perturbations 
may be calculated by harmonic analysis. The breaking of the 5*0 (4) isometry at the tip 
of the KS throat has been studied in [211] (see also [44!|), where it was interpreted in the 
language of the AdS/CFT correspondence. From the point of view of the field theory dual 
to the warped background, the gluing of the compact CY to the throat may be described 
by the addition of symmetry breaking operators in the UV. The RG flow of these operators 
corresponds to the radial profile of the perturbations of the 10-dimensional background. 
On an AdS^ x background, a perturbation eiy) Sfij^fl^) is decomposed into harmonics 
5fi^\Vl^) on X^ each with its own profile in the radial direction 

e{y) 5^(fi5) - (^5) , (12) 

where we have neglected the sub-leading contributions with respect to y. The radial profile 
follows from the linearized radial mode equation after dimensional reduction on . The 
parameters Ul thus depend on the mode under consideration and its quantum numbers L 
on X^. The leading perturbation in the throat arises from the mode with the lowest value 
of aL- 



Ref. 2l| classified the perturbations, breaking the 5*0(4) isometry but preserving 4- 
dimensional Lorentz invariance and supersymmetry, which may be turned on a KS throat 



approximated by AdS^ x T^'^. The leading perturbation found in 2l|] corresponds to a 
perturbation e'"'^^^ ^ 6 f-^ of the metric with ai = v^28 — 4 ~ 1.29. It involves a tensor 
mode 6f^^\Q5) on T^'^ with quantum numbers L = {li = l,l2 = = 0), where {li,l2,r) 
are the principal quantum numbers associated to the SU{2)i x SU{2)2 x f/(l) isometry of 



T^'^ (see 22] for the KK spectroscopy of the type IIB supergravity fields on AdS^ x T^'^). 
The isometry breaking perturbations involving other harmonics are exponentially suppressed 
in the throat compared to the leading one. Note also that, for a perturbation with given 
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radial profile e~"^^/^, the principal quantum numbers are fixed but the other ones (those 
which are not involved in the eigenvalues of the Laplacian on X^) may vary since they do 
not appear in the radial mode equation. 

Finally, Ref. j21| estimated the perturbation of the warp factor e^"^ at the tip of the KS 
throat induced by the perturbation 6fij. The leading correction was found to be of the 
order 6A(yt) ~ Q-^vt/^ with the same coefficient a ~ 1.29. This mode is absent on the 
AdS^ X T^'^ background but should correspond to a linear perturbation around the full KS 
solution. Again, we expect different harmonics in the angular directions to have different 
exponential suppressions at the tip of the throat. Because these modes arise from deviations 
of the KS solution from AdS^ x T^'^, their radial profile in the throat may differ from e~'^yl^. 
However, what will be important for us is their suppression at the tip, bA[y^) ~ ^-(^yt/R^ 
and it will be convenient to keep the profile e~"^/'^ in the throat to model this suppression. 
The perturbation of the metric in the angular directions and the perturbation of the warp 
factor will play similar roles in the following. 



C. Parameters of the Warped Geometry 

Here we discuss the basic parameters of the warped compactification, which will enter into 
the formulas for the cosmological abundance and decay time of KK modes. The cosmological 
constraints from the dangerous KK relics will be given in terms of these parameters. 

One of the essential parameters will be the leading exponent ai (for the KK modes with 
quantum numbers L) for the isometry breaking, entering in ( ITIl) . ( fT2l) . In addition to the 
parameters a, and to the string length ^fa' and coupling Qs-, our treatment will involve 
essentially three other parameters: the (AdS) radius of curvature i? in ([8]), the volume of 
the whole 6-dimensional internal space Vg, and the warp factor at the bottom of the throat 
Q-yt/R_ w/"g now recall useful relations between these parameters. 

The volume Vq is dominated by the "big" compact CY manifold, and it is always greater 
than the contribution of the throat alone, V^^^ > R. For instance, for the KS throat this 
implies Vq > j^^^^^^d^yVd ^ Vo\{T^'^)R%/2 = M^R%/27. We will also take R > ^ ivl 
order to work in the supergravity approximation. In KS, R > R- > \/a' may be achieved 
by a suitable choice of fluxes, see Eq. ([7]). 
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The string length and Vg are related through the value of the 4- dimensional Planck mass 

^ wrk^' ■ 

We will take Qs = Q.l in the estimations below. 

The final parameter is the value of the warp factor at the bottom of the throat, e~^*/'^. 
We will distinguish the short throat and the long throat, when this parameter is relatively 
big or relatively small, correspondingly. The short throat is associated with the warped 
brane infiation proposed in 6], and is often called the infiationary throat. In the original 
model, the normalisation of the adiabatic density perturbations generated during inflation 
reads as 

5h ^ 0.4Ar5/6 3 g-4y,/3iJ ^ 2 X 10^^ (14) 

\Mpi/ 

at a number of e-foldings Ne ~ 60. Here T3 = [{2ny^gsa''^] ^ is the tension of the D3, D2> 



branes responsible for inflation. The origina 



generalized by several models (e.g. , llO, lUl]) for which the parameters may be different. 
Eq. f[T^ will serve as a convenient reference value for a short throat. 

The long throat is associated with the hierarchy between the string scale and the TeV 
(electroweak) scale a la Randall-Sundrum js], if the Standard Model degrees of freedom are 
located on a brane around the tip of throat. In this case the warping is signiflcant 

(15) 

where Mg = 1/ v^- However, it does not mean that inflation cannot be associated with the 
long throat. The original KKLMMT throat slow roll inflation suffers from the ?7-problem 
due to the conformal coupling of the inflaton. It turns out that one still can obtain sufficient 
fast roll inflation with the conformal coupling, but only in the long throat 8|. For this 
case, the estimation f[T^ for the amplitude of the inflaton cosmological fluctuations is not 
applicable. 

Both Eqs. ( |T5i) and ( |T^ determine the value of the following speciflc combination of the 
parameters 

I ey^/^ . (16) 



KKLMMT inflationary construction has been 



'a' 

We will thus use this combination in the calculations of the KK modes abundance and decay 
times. 
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For the short throat (fT^ . using Eq. (fT3l) with Qs = 0.1, Eq. (fT6|) gives 

A/" 5 X 10^ ( short throat) . (17) 

In this case, significant warping (e^'/^ ^ 1) requires at least Vq^^/^/o' <^ 80, and smaller 
values for R/y/a' and V^^^/R. For example, Vg^'^^ ^ 5 leads to e^^^^ ^ 4 x 10^. 
For the long throat, using Eqs. ( fT3ll and ( fTSl) . we have 

AT ^ 10^^ (long throat) . (18) 



III. SPIN 2 KALUZA-KLEIN MODES 

Let us recall the energy cascade at the end of the brane-antibrane inflation [l^. Brane- 
antibrane annihilation through tachyon decay results in the production of massive closed 
string loops, which further decay into KK modes, interacting with each others, and with 
SM particles. Massive KK modes are localized at the bottom of the throat, while massless 
four-dimensional gravitons are not. There are various KK modes of the inner manifold 



with the throat (s) with fluxes, e.g. [23|, |2J]. We will be interested in the KK modes with 
angular degrees of freedom associated with inner dimensions. In the KK sector with angular 
momentum, bearing in mind cosmological applications, we shall be primarily concerned with 
the lightest modes localized at the bottom of the throat. Unfortunately, the mass spectrum 
of KK modes in the KS geometry is not known exactly but only up to the order of magnitude 
estimation m ~ e~^^f^/R. 

At this point the study can go into different directions, depending on the choice of the 
prototype KK relics. Ref. considers a scalar KK mode, which is the lightest one on the 



AdS^ X T^'^ background ^. However, to make this mode unstable, [27| introduced additional 
SUSY breaking operator, and the final decay products involve massless gravi-vector modes 
^. There can be other subtleties: next-to-the-lightest modes may be forbidden to decay 
into the lightest modes, and their decay time directly into SM brane degrees of freedom 



Such a mode is not necessarily the hghtest one after compactification (for this mode, specific boundary 



conditions have to be imposed on a sHce of AdS, see 25|, |26[) and after moduH stabihzation 



^ In the third version of [23| , the decay channel has been modified and the final decay products involve only 
massless gravitons. 
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may be different. We will encounter these situations in our calculations below. Therefore, a 
comprehensive investigation shall include different possibilities. 

In this paper, as a representative set of these KK modes, for economy and definiteness, 
we focus on the spin-2 fields in four-dimensions, and their decay into excitations of the 
SM brane located in the same throat. The dynamics of the spin-2 KK modes follows from 
the ten-dimensional linearized Einstein equation. Their wave equation depends only on 
the 10-dimensional background metric and is decoupled from the fluctuations of the matter 
sources. Our purpose is to calculate the abundance and lifetime of angular KK modes 
and make connection between cosmological constraints and parameters of the model. The 
methodology we develop in this paper can be extended for different KK modes and different 
phenomenological settings of their decays. 

The spin-2 fields correspond to the symmetric transverse-traceless perturbation /i^i, 

g^.u = rit^y + h^,{x^, with r^^'^^A V = V = • (19) 

of the 10-dimensional metric ([2]). 

The equations of motion for the KK modes fll9p in the background geometry ([2]) are de- 
rived in the Appendix. The 4-dimensional Poincare invariance allows to solve by separation 
of variables 

h,,{x\ y^) = J2 ^miy") 7i:^n^') , (20) 

m 

where 7^™'*(x) are the purely 4-dimensional spin-2 fields of mass m, satisfying 

□(4)7i^^ = m^i™) (21) 

with 0(4) = Tj^^d^dy. The wave equation for the radial profile of the modes in the internal 
dimension may be written 

-Ve(e^^V^$,n) =m2e2^$^, (22) 
where Vc is the covariant derivative associated to the metric in the internal space, cjab in 

The spectrum and profile of the modes are mainly determined by the region where the 
warp factor e^^ varies exponentially. We will first consider the KK modes in the local 
isometric background ([8]). When isometrics are broken according to (ITT]) , the metric per- 
turbations are suppressed, e{%)) ~ er'^'^l^ <^ 1, in the main part of the throat. We will solve 
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the KK mode equation around this non-isometric background at hnear order in e, by using 
the standard techniques of quantum mechanics perturbation theory (sub-section IIII B| ) . 



A. KK Modes of Isometric Throat 

To study the spin-2 KK modes around the local geometry ([8]), we may use the results of 
the Appendix with e^'^^y"^ = e~^^/'^ and Qabdy'^dy^ = dy"^ + fij dO'^dOK In this geometry, 
the mode equation fl2^ reads 



d_ 

dy 



-Ay/R 



d^ 

dy 



-iy/R 



i?2 



m 



2 e-2s//« $ 



(23) 



where Af denotes the Laplace operator associated with the metric /^(^^s) of the angular 
coordinates. 

The background symmetries lead to the separation of variables 



$m(/)=^nL(l/)Qf(^^5 



with 



AfQi' = -F'\L)Q 



M 
L ; 



(24) 



(25) 



where L denotes collectively the quantum numbers involved in the eigenvalues F'^{L) of Aj, 
and M the other ones. For instance, if fij is the metric of the 2-sphere, are the usual 
spherical harmonics Y^^ and F"^ = L{L + 1), where L are positive integers and M are integers 
with \M\ < L. For a Klebanov-Strassler throat, Qf^ would be the harmonics associated to 
the angular dimensions of the deformed conifold. 

Substituting fl24|) in Eq. (123|) and using (125!) gives the radial wave equation 



d 
dy 



-Ay/R 



dip. 



riL 



dy 



-Ay/R 



+ 



F^{L)i)nL = ml^e 



-2y/R 



(26) 



whose eigenfunctions ipni and eigenvalues m^^ will depend on the quantum number(s) L 
(through F'^{L)) and on a radial quantum number n. For the massive modes (m^ 7^ 0), the 
solution for the radial wave-function is 

^nL(2/) = iVnLe'^/^ [j.(m„ii?e^/^)-5„in(m„ii?e^/^)] , u=^A + F^{L) , (27) 

where J^, and Y^, denote the Bessel functions of order u of the first and second kind. The 
constant BnL and the mass spectrum ninL are determined by the boundary conditions, while 
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the constant NnL is fixed by the normahzation condition. In general, KK modes are strongly 
localized at the bottom of the throat, y = yt, and their masses are of the order m ~ ^ 



Boundary conditions may be obtained by matching the solution ( |27I) to the solution for 
the modes at the base (at y = 0) and at the tip of the throat (at y = yt). The effect of the 
CY on the boundary condition is difficult to model since we don't know the explicit metric 
of the gluing region. Here we will focus on the modes in the interval < y < yt and for 
simplicity impose Neumann boundary conditions at both ends. To the order of magnitude 
we will be eventually interested in, we don't expect the results to be very sensitive to the 
particular form of the boundary conditions. The conditions ipnLi^) = i^'niiyt) = translate 
into two conditions for the Bni 

{v - 2) J^(m„Li?) - rUnLR Ju-i{mnLR) 



{v - 2) Y^{mnLR) - m„LRyu-i{rnR) 
{v - 2) J,(m„ii?e^*/«) - m„i/2e^*/« J,_i(m„Li?e^*/^) ^^s) 



(z/ - 2) Y,{m^LRey^/R) - m„i/2e^*/« n_i(m„ii?e?''/«) ' 

In the following, we will be interested in the lightest among massive KK modes, for 
which mR <^ 1 . In this case, the first equation above implies Bni << 1. More precisely, 
BnL oc {mRY for z/ = 2 and BnL oc {mRY^ for z/ > 2. The second equation in (l28l) then 
reduces approximately to 

(2 - V) .UU) + U Ju-liCnL) = , (29) 

which determines the spectrum of the KK eignemasses. We have defined ^nL according to 

m„i = ^ e-y^/^ , (30) 



where n = 1,2,3,.... counts the successive roots of Eq. fl29l) for a given z/. The mass of 
the KK modes increases with n and with u, i.e. with F'^{L) in (|25l) . For S,nL ^ 1 we find 
an analytical approximated solution (obtained by expanding the Bessel functions for large 
argument) 



Some KK modes of mass m ^ Vq may be localized instead in the bulk of the CY, see e.g. [38]. We will 
not consider such modes here because for sufficient warping they are much heavier: 3> er^*^! ^' j R. 
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When the throat is approximated by a patch of AdS^ x T ' , we can use the eigenvalues 



33 



3^. We have 



of the Laplacian for the scalar harmonics fl25|) on T ' , see e.g 

F^{L) = 6 (^hih + 1) + hih + 1) - ' (32) 

where L = {li, I2, r) are the Sf/ (2)i x Sf/ (2)2 x f/(l) principal quantum numbers, with li and 
Z2 both integers or both half-integers, and with r/2 G {— and r/2 G {— ^2, ^2}- 
In this case we can calculate the spin 2 mass spectrum of the angular KK modes. For 
illustration, we list in Table I the lightest massive states when the throat is approximated 
by a patch of AdS^ x T^'^ with Neumann boundary conditions for the radial wave functions. 
They are tabulated according to their radial quantum number n and their angular quantum 
numbers li,l2,r- 



n 


h 


I2 


r 


F2(L) 




1 














3.83 


1 


1/2 


1/2 


1 


33/4 


5.45 


1 


1 








12 


5.98 


1 





1 





12 


5.98 


2 














7.02 


1 


1 


1 


2 


21 


7.04 


1 


1 


1 





24 


7.35 


1 


3/2 


1/2 


1 


105/4 


7.57 


1 


1/2 


3/2 


1 


105/4 


7.57 


1 


3/2 


3/2 


3 


153/4 


8.63 


2 


1/2 


1/2 


1 


33/4 


8.92 



TABLE I: Lightest massive spin 2 KK modes of the throat approximated by a patch of AdS^ x T^'^ 
with Neumann boundary conditions for the radial wave functions. 



Now we compute the normalisation constant NnL- The normalisation condition for each 
'y'iiu\x) in (1201) to represent a canonically normalized spin-2 field from the 4-dimensional 
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point of view is given in Eq. (11071) . We normalize the wavefunctions in the angular directions 
according to 

J d^n, Qf (fis) Qf>'m = hu hiw , (33) 

where dfl^ = dO^ ... dO^ \ff is the element of solid angle. Eq. fll07p (here for D = 10) then 
gives 

^ / e-'y/^^lj^M^n'Mdy = 6^^, (34) 



4 

for the radial profile of the modes, where Mf^ = j2n)^gW^- range of masses we are 

interested in (e~^*/^ < mR <^ 1), this integral is dominated by the upper limit, where we 
can neglect the term in Bniyu- This gives 



(35) 



see e.g. 



361]. It is more convenient to rewrite this in terms of the 4-dimensional Planck 



mass, by using Eq. fllOSp : 



where we have also neglected the last term in parenthesis (which is approximately equal to 
one). 

For the modes with angular momentum, the second term in (127|) is comparable to the 
first one at ?/ = and becomes quickly negligible as y increases. Their normalized wave 
function is thus well approximated by 



Finally, Eq. fl26l) has a solution with m^^ = if F'^{L) = 0. This mode, with zero mass 
and no angular momentum, is the usual 4-D graviton. Its wave function is constant in the 
extra dimensions and normalized to 

$o = ^ooQS = t^, (38) 
Mpi 

see ( 11091) (for this specific mode, the the leading contribution to the normalisation constant 
comes from the bulk of the CY) . 
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B. KK Modes of the Throat with Isometry Breaking Perturbations 



We now study the spin 2 KK modes in the background geometry (fTTj) . where the isometry 
breaking perturbations are suppressed by e{y) = e""^/"^ in the throat. The corresponding 
mode equation is given by Eq. (12^ with e'^^^''^"^ = e'^^/^ [1 + e(|/) ^(^5)] and Qabdy^'dy'^ = 
dy"^ + [fiji^Q^) + e{y) Sfij^Q^)]. The function e{y) is much smaller than one in the throat 
itself and we can treat it as a perturbation. Up to first order with respect to e, the mode 
equation may be written as 

[Ho + V] <l>™ = m2e-2^/^$™ , (39) 



where 



° dy 



9 ^ _ (40) 



dy 

is the unperturbed operator involved in the LHS of the mode equation (123|) for the isometric 
background, while 



V = ewHo - e-^^/«^ (2w + -f^Sh!] # 
" dy \ 2-' -""y dy 



is the perturbation operator encoding the effects of isometry breaking 



(41) 



To solve for Eq. (!39l) . we may follow the standard procedure of perturbation theory, 
by expanding the wave function on the complete set ipniiy) Ql {^5) of solutions of 
the unperturbed problem of the previous sub-section, and determine the coefficients of the 
expansion by using the orthonormal conditions (!33| [34|) . We then define the matrix elements 

Kl^l' = / dydn,ij^L{y) (l^s) ^ ^n'L'(2/) Qui^n) (42) 

of the perturbation operator V in the basis {ipniQ^}- The eigenvalues at zeroth order in e, 
i.e. the KK masses squared m^^ in the isometric background, will typically be degenerate 
since in particular they do not depend on the quantum number(s) M of the harmonics 
QfJlQ^). In this case, the correct wave functions at zeroth order are the linear combinations 
which diagonalize the matrix elements of degenerate states (see e.g. 37|]): V^^'y = for 
T^nL = rrin'L' and {n,L,M) 7^ {n',L',M'). We will keep the notation iPulQ'l these 
correct wave functions at zeroth order. 
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The KK modes' masses up to first order in e are then given by 

rnlLM = ^Il + \ Ml,R' Kl^l , (43) 

where irinL are the masses at zeroth order considered in the previous sub-section. In general, 
the breaking of the isometries removes the degeneracy. We will be specially interested in 
the wave functions of the KK modes when the isometries are broken. Up to first order in e, 
they are given by ^ 

1 A/f® yM'M 

n',L'.M' ^ ^ 

(45) 

where the ellipses stand for similar contributions due to the matrix elements between de- 

n 

generate states (see [37l]). 

Of course, there is still a zero mode 4-D graviton, with zero mass and constant wave 
function, in the absence of isometries in the internal space. Indeed for this mode, the 
corrections to the zeroth order mass and wave function vanish: V^[qq = because V'ooQo is 
constant. 

Consider now the matrix elements ^^^'/^^ between two massive states. The contribution 

of the first term in (|4TI) to (l42l) may be written 

ryt r 
Ka''nL = <Ly^ dye-'y/^e{y)^nL{y)^n'L'{y) J dn,w{n,)Q^{n,)Q^:{n,) + .... (46) 

The other terms in (HT!) lead to similar contributions. The integrals over the angles deter- 
mine the transitions allowed (i.e. the non- vanishing matrix elements), while the y-integral 
determine their amplitude. 

In the following, we will need the amplitude of the matrix element between two light 
massive states {mnL ~ ^n'L' ~ e^^'^^/i? ^ 0) with non- zero angular momentum (L, L' ^ 0). 

^ The matrix ()42p is not symmetric, but one may derive the relation 

V^f^^^L ~ y:^k' - - <l) J dy dn, e-^y'^ e [w + Sf,^ Qf i^wL' Q^' • (44) 

From this it can be shown that the wave functions (|^5l) satisfy the orthonormal conditions (|107p at first 
order in e, as they should. This fixes the constant in 0{e) in the first term of (pS]). but its expression will 
not be important for what follows. 
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Using the expression flTTI) for the tpni and with the change of variables u = rjinLR^^^^^, the 
y- integral in (H6il for a perturbation e{y) = e~"^/^ gives 







where we have used ( l30l) and ^ni ~ C'(l) for the hghtest modes. For m ^ 0, Jy{u) oc m'^, so 
the u-integral does not depend on its lower limit for a < 2 + u + u' }^ The integral is then 
mainly accumulated around its upper limit, leading to the last equality in P7|) . 

Therefore, for the non- vanishing matrix elements V^^!^^ between two light massive modes 
with angular momentum, the component of the first order wave function in (H5l) along 
the zeroth order wave function ipn'uQ^' is given by 



1 A/fS R5 yM'M 



PI 



where we have used fllOSp and (!30|) . and keep the contribution of one of the terms in the 
sum (1451) . Note that the first order correction to the wave function is smaller than the wave 
function at zeroth order by a factor of t{yt) at any position in the radial coordinate y. 

IV. INTERACTIONS OF KK MODES 

In this Section, we will study interactions of the KK modes. We shall include 3-legs and 4- 
legs self-interactions of the KK modes. Most important for cosmology will be interactions of 
KK modes with the Standard Model particles. This will depend on the SM phenomenology 
in the string theory setting. Following [l^, we will consider a simple model of KK modes 
interacting with the SM located on a 3-brane around the tip of the same throat. Notice 



that this model differs from what was considered in 



27l |. The cases where the Standard 



For a > 2 + V + v' ^ the u-integral diverges for u — > and is then dominated by its lower limit, leading to 

vt , ^ ,^ -(2+u+u')vt/R y„ 

dye~^'+"^y^/^i^Mi^n'L'{y)^ (48) 

This case may be covered by replacing a with Min(a , 2 + z/ + z/') in Eq. ([^5]) . 
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Model degrees of freedom are located in the bulk of the CY or in another throat are further 
constrained and will be discussed in the last, conclusion Section. 

The decay times of angular KK modes due to isometry breaking are much longer than 
the time scales for their annihilations. It is convenient to call them as the "late" decays due 
to isometry breaking versus "early" decays and annihilations which are taking place without 
isometry breaking. Therefore, we consider first the "early" channels of KK interactions, for 
which the isometry breaking of the throat can be neglected. At this level we will deal with 
KK modes with conserved angular momenta (of inner isometrics). 

We first compute the coupling constants for KK modes self-interactions in sub-section 
IIV A[ In sub-section IIVB[ we then compute the coupling constants for the KK modes 
interactions with the SM at the brane, again, with no isometry breaking from the CY. In 
both cases, we pay special attention to the selection rules to be satisfied. This will allow us 
to identify the long-living angular KK modes, which could pose cosmological problems. We 
then study in sub-section IIV CI the "late" decays of the long-lived KK modes induced by the 
isometry breaking perturbations. 

A. Self-interactions of KK Modes without Isometry Breaking Perturbations 

For the decays and annihilations of KK modes into themselves, the effective four- 
dimensional Lagrangian and coupling constants are obtained by substituting the pertur- 
bations (fT9|) into the bulk action, expanding to third and fourth order in h^^, and integrat- 
ing over the internal coordinates. This is done in details in the Appendix for the general 
background ([2]), up to third order in h^i,. 

Here we consider KK modes self-interactions which are already present in the throat 
geometry ([8]) without isometry breaking, i.e. at zeroth order w.r.t. the isometry breaking 
parameter e{y). 

We start with the decay of a KK mode of mass m and quantum numbers (n, L, M) 
into two lighter KK modes, of masses m' and m", and quantum numbers (n', L', M') and 
{n'\ L", M"). The corresponding terms in the effective action have the form 




(50) 
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where dQ^ = dO^ ... dO^ \fj is the element of sohd angle in the angular dimensions X5. The 
expression above may be obtained by substituting ( !24|) into ( lllOp for the background (JHl). 
Here M\q is the 10- dimensional Planck mass, which is related to the 4-dimensional Planck 
mass Mpi and the total volume of the internal space Vg through Eq. (11081) . With our 
choice of normalisation for the wave function in the internal space, the 7^™^ are canonically 
normalized spin 2 fields, and Askk is the corresponding coupling constant. Eq. (!50|) allows 
us to calculate Askk- 

Let us consider the selection rules that follows from Eq. ( !50l) . The conservation of the 
4-dimensional energy and momentum in the center of mass frame reads as 



+ k2 + Vm"2 + k2 , (51) 



which requires in particular m > m' + m". 

The integral over dfis imposes angular momentum selection rules. By contrast, the 
warped background has no translation invariance in the radial direction so the y-integral 
does not impose the conservation of the total KK modes mass: m may be different from 
m' + m", contrary to the case of extra dimensions with translation invariance. However, if 
one of the decay product is the massless graviton, m" = L" = M" = 0, its wave function is 
constant in the internal space and drops out of the integrals over the inner dimensions. With 
( l33l) . the y-integral then reduces to the orthogonality conditions which implies m = m'. 
Furthermore, decay of KK modes is cascading into lighter and lighter KK states. However, 
if one of the decay product is the zero- mode graviton, m" = 0, then m = m' is required, 
as we just saw above and is shown more generally in the Appendix (see Eqs. (IllOp and 
(lllip ). Therefore, a massive KK mode (jn^ 0) cannot decay into two massless gravitons 



m' = m" = (this also was noticed in [l5l])- It is also obvious that a massive KK mode 
cannot decay into one massless graviton and another massive mode since the condition 
m = m' leaves no phase space for this process. Therefore, a massive KK mode cannot decay 
by 3-legs interactions into one or two zero-mode graviton(s). As shown in the Appendix, 
this is independent of any isometrics in the internal space. In addition, suppose mo is the 
mass of the lightest massive KK mode (mo 7^ 0). Evidently, only the modes with m > 2mo 
can decay into two other massive modes. We thus conclude that the modes with m < 2mo 
cannot decay by 3-legs interactions into any other spin 2 modes (massive or massless). 
Finally, in the presence of isometrics, the conservation of the corresponding angular mo- 
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mentum imposes further selection rules. In particular, a mode with non-zero angular momen- 
tum (L 7^ 0) cannot decay into modes with non-zero angular momentum only (L' = L" = 0), 
which implies that the lightest modes with non-zero angular momentum cannot decay into 
any other KK mode. 

For the allowed interactions between 3 massive KK modes, the effective 4-dimensional 
coupling constant is obtained by substituting Eq. (1371) into the l.h.s of Eq. (150|) . and per- 
forming the integrals over the internal coordinates. The ^/-integral is accumulated mostly 
around y = yt, leading to 

where we have neglected numerical factors of order unity and used ( llOSp to express M^q in 
terms of Mp\ and Vq. 

The coupling constants for KK modes 4-legs interactions are obtained in the same way as 

we did it above for the 3-legs interactions. Expansion of the bulk action to the fourth order 

w.r.t. h^^ gives the effective coupling constant for interactions between four KK modes 

ryt 

A4KK ~ MIqR^ I e"^^/^ i)nL i^n'L' i^n"L" i^n"'L"' dy . (53) 

Jo 

For four massive modes, the integral is again dominated by the contribution around its 
upper limit, leading to 

A4KK-^r^. (54) 



It is worth noticing that this is about the square of fl52|) . A4KK ~ ^Ikk- 



B. KK-SM particles Interactions Without Isometry Breaking Perturbations 

To obtain the coupling constants for the decays and annihilations of KK modes into 
SM brane degrees of freedom, we consider 4-dimensional fields located on a probe brane at 
y'^ = yl in the geometry ([2]). One can start with the DBI action for the probe brane and use 
the decomposition of the DBI action for the perturbative field fluctuations. For simplicity, 
for a light scalar four- dimensional brane degree of freedom if, we have 

j d^'x v/^G^Gj:;,) d.Hd^H 5(^)(/ - yl) = j d'x -g{yl) g^^'iyl) d.Hd^H , (55) 

where G^^^f = e^^'^'^b) g^^[yl) is the induced metric on the brane located at y^ = y^. To zeroth 
order in hfj,^, the kinetic term is canonical for the normalized scalar field H = e^^^^) H, that 
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we have substituted in the r.h.s. of Eq. 0551] . The interaction term h^^ d^H d^H is then 
obtained by substituting g^^ = rj^i, + h^^, where hfj_^ are the spin 2 KK modes f|T9l) . Later 
on, we will consider the branching ratios for the KK modes decay into different Standard 
Model fields. For the moment we are only interested in the overall coupling constant, which 
can be estimated for the simplest case of a light scalar field. 

Consider the decay of a KK mode with quantum numbers {n, L, M) in ([21]) into two 
lighter 3-brane scalar degrees of freedom, KK — b b, following from (1551) . It involves the 
term 

i^uM j d'^x^i:^^ d,Hd,H = AxKbb j t^'x7f^) d,H d,H , (56) 

where fib and yb denote the position of the brane in the angular and radial dimensions 
respectively. 

At first glance, the coupling constant Akk6 6 somehow depends on the value of the angular 
eigenmode of the KK wave function at the position of the brane. While the brane will 
break some, but not all, the isometries of the bulk, the interaction (!56|) should respect the 
conservation of angular momentum associated with the remaining isometries. 

To illustrate that the interaction in the form (|56ll indeed conserves quantum numbers 
associated to the isometries unbroken by the brane, we consider an example where the bulk 
involves an exact 2-sphere 5^. The setup is shown in Fig. [TJ The angular dependence of the 
KK mode wavefunction on S"^ is given by the usual spherical harmonics, Q^^ oc ¥^'^{6, 0). A 
3-brane breaks the 5*0(3) isometry of the S*^ into 50(2). Choosing, for instance, the 2;-axis 
to point in the direction of the position of the brane on the S"^ (here the north pole, see Fig. 

9b = 0, the coupling constant in ( l56l) is proportional to Akk6 6 oc Y^^iOb = 0) oc Pl^{0) 
(where P/^ are the associated Legendre polynomials [36|). Notice that this vanishes for all 
M 7^ 0, see Fig. [1] (where M is the projection of the angular momentum along the 2;-axis). 
This means that only the modes with M = can decay. The conservation of the quantum 
number M associated to the isometries unbroken by the brane forbids the modes with M ^ 
to decay. Indeed, for these modes, the angular momentum vector has a component along 
the z-axis, and rotations around this axis are left invariant by the brane. Analogously, a 
3-brane on the 3-sphere at the tip of a Klebanov-Strassler throat breaks the 5*0(4) isometry 
into 5*0(3), and KK modes with non-zero angular momentum associated with the remaining 
isometries cannot decay. 
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SM D3 branes. 




Zeros and signs of sperical harmonics at S 

FIG. 1: Zeros and signs on the spherical harmonics on 5^. We compare the M = with the M 7^ 
case. The spherical harmonics with M 7^ vanish at the two poles. As a consequence, KK modes 
with nonvanishing angular momentum along the directions whose isometries are left unbroken by 
the brane are not directly coupled to brane fields (see the main text for details). 

For the allowed interactions, the coupling constant AxKbb in ( |56|) is obtained by evaluating 
the wave function (l37j) at the position of the brane, y = yb- For the decay of a massive KK 
mode {rrinL 7^ 0), this gives 



X 



KKbb 



(57) 



J Mpi J^imRey^l^) ' 

Comparing with Eq. (!52|) . we see that the 3- legs coupling of a KK mode with brane degrees 
of freedom located at the tip of the throat, yb = Vt, is of the same order as the 3-legs coupling 
of the KK modes between themselves. We will denote this coupling constant by A3 in the 
following 

AsKK ~ AKKbb ~ ( 1 = A3 . (58) 

Interactions involving two KK modes and two brane degrees of freedom are obtained 
by expanding (|55|) to the second order in h^y. The corresponding coupling constant is 
obtained in the same way as in (!56l) . Again, it vanishes for KK modes with non-zero angular 
momentum associated with the isometries unbroken by the brane. For two other, massive 
KK modes, it is given by 

y\ ^2{2y,-yt)/R J ^[mRe^^ > ^) J y> {m' Rcy^ > ^) 



A 



KKKKb 



b ~ 'ipnliybj'ipn'L'iyb) 



(59) 



which again reduces to (IMI) when the brane degrees of freedom are located at the tip of the 
throat, yb = yt- We will denote this couphng constant by A4, and it is related to A3 in ( 158]) 
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according to 



A4KK ~ AKKKKbb ~ A3 = A4 



(60) 



C. Late Decays of KK Modes due to Isometry Breaking Perturbations 

In the isometric throat, the KK modes with non-zero angular momentum associated with 
the isometries unbroken by the brane (the modes with M 7^ in the 5*^ example) cannot 
decay into the brane degrees of freedom. Among these, the lightest ones cannot decay 
either into any other KK mode. They are therefore stable if the isometries are not broken, 
and become long-lived when the isometries are slightly broken. This is illustrated by the 
diagrams of Fig. [2l 

We now study how these long lived modes may decay in the background ffTTj) . where the 
isometry breaking perturbations from the CY have radial profile e{y) <^ 1 in the throat. We 
proceed as in the previous sub-sections, but now we consider the effective action at first order 
in e. The corresponding wave functions for the KK modes were calculated in sub-section 

iniB 




k 



FIG. 2: Left panel: for the isometric throat, the decay of angular KK modes localized around the 
tip of the throat is forbidden by higher-dimensional angular momentum conservation. Right panel: 
The isometry breaking leads to a new interaction, which corresponds to the decay of an angular 
KK mode into SM particles mediated by a background isometry breaking KK mode localized in 
the compact CY. 

Consider first the interactions between 3 KK modes described by flllOp . The correspond- 
ing 4-dimensional coupling constant at first order in e includes terms involving the first order 
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correction to e^^ and ^/g, and 3 zeroth order KK modes wave functions 
A?L = / dye-^y"'e{y)ilJnL'ipn'L''ipn"L" j dn, (^ + }-f^'5h^ Ql Qu Qu' + - ■ (61) 

In addition, there are terms involving the unperturbed background metric and the first order 
e correction to the wave-functions, not shown in fl6T]) . The complete expression is given by 
Eqs. dnOD and ffTTTD . 

In order to evaluate Ag^^j^ it is sufficient to consider the first term in ( 16T|) . There is, 
however, an exception. While the first term in fl6Tl) formally does not forbid to have massless 
m = KK mode in that three-legs interaction, the total amplitude in this case is zero. 
The reason is that the trilinear combination of the full KK wave functions in Ag^^^j^, when 
one of them is the constant zero-mode wave function, is effectively reduced to a bilinear 
combination. This bilinear combination is constrained by the orthonormality conditions 
of the full wave functions, which forbids the three-legs interactions involving one or two 
zero-mode (s). See the appendix for details. 

In general, the isometry breaking perturbations are decomposed into harmonics w'''-'\VL^) 
and flp{^l^) with different quantum numbers £, each with its own profile in the radial 
direction, see section Hi Bl and Eq. ( fT2l) . Consider the selection rules resulting from a given 
harmonic f^p with radial profile e{y) = e'"^"^^^ in ( |6T1) The function f^^^fl^ in the 
integral over the angles breaks the conservation of angular momentum, which may facilitate 
the decay cascade of the angular KK modes. However, for the lightest of them, the relevant 
process is their decay into modes with zero angular momentum (L' = L" = 0), which 
requires / dn f^SfipQfl ^ 0. The perturbation f^6fiP having given principal quantum 
numbers C, it allows only the modes with the same angular momentum, L = C, to decay 
into two modes with zero angular momentum. The lightest modes with angular momentum 
would then decay through the perturbations carrying their own quantum numbers. These 
perturbations may be exponentially suppressed compared to the leading one, resulting in 
exponentially longer lifetimes. Furthermore, as in sub-section IIV At the decay of the modes 



The leading perturbation Sf^f'' (the one with the lowest ac) on AdS^ x T^^^ identified in [2l| corresponds 



to a traceless tensor mode on T^'^, for which f'^^Sf^f' vanishes. In this case, the contribution of the 
leading perturbation to Ag'^^j^ comes from the perturbed wave- function (through the last term in (j4ip ) or 
from the perturbation of the warp factor, leading to similar results. Here we consider the contribution of 
another harmonic for illustration and keep ac arbitrary. 
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with m < 2 mo is still kinematically forbidden, since this does not depend on the isometries 
of the internal space (see the Appendix). As illustrated in Table 1, there may be several of 
these modes. Because of these restrictive selection rules, the leading decay channels for the 
long-lived modes will be into brane degrees of freedom. 

For the decay channels which are made possible due to a given isometry breaking pertur- 
bation with radial profile e~"^/^ in the throat, the coupling constant Askk is calculated by 
evaluating the ^/-integral in fl6T]) as we did in fH7|) . For 3 light massive KK modes, this gives 

^--(1) (62) 

We consider now the decay of a KK mode into 2 brane degrees of freedom at first order 
in e. We are interested in the long-lived KK modes, i.e. the lightest modes with non zero 
angular momentum associated to the isometries unbroken by the brane. As in fl56l) . the 
corresponding coupling constant is just given by the wave function evaluated at the position 
of the brane in the internal space. This vanishes at zeroth order for the modes we are 
interested in, so we have to consider the first order corrections to the wave functions $^ 



given m 

The decay occurs through the combined effect of isometry breaking due to the CY and 
isometry breaking due to the D?) (or Z)3) brane. The resulting selection rules are less 
restrictive but more involved than for the decay into two KK modes with zero angular 
momentum. As an illustration, we return to the example of a 3-brane located at = on 
a 2-sphere {Q^ = Y^'^). The modes which could not decay into brane degrees of freedom 
in the isometric case were the ones with non-vanishing quantum number M. In the series 
f H5l) . all the terms involving Y^f' with M' ^ still vanish at the position of the brane. The 
4-dimensional coupling constant is thus given by 

= \ E ^^"^"S Yl,{9, = 0) + ... . (63) 

n .L 

This coupling constant is nonvanishing if the wave function acquires a non-zero 

component along some harmonic Y^, due to the perturbation, i.e. if there exists some L' 
such that V^^li^i^ 7^ for M 7^ 0. Using (H^ . the contribution of the first term in fHTi) to 
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this condition may be written 

j w Yl^ Yl, ^0 for M . (64) 

Consider the contribution of the leading perturbation (the one which is the less suppressed 
at the tip of the throat) w^'-'^ with given principal quantum numbers C In our 2-sphere 
example, we would have w^'-''> = J2m=-c^m'^c^ ■ this case, the f22-integral in (IMI) is 
nonvanishing if there exists an L' such that \C — L\ < L' < C + L with C + L' + L even, and 
if there exists an Ad such that = M. The first condition can always be satisfied, while 
the second one requires L < C. Thus in this case the KK modes can decay into SM (brane) 
degrees of freedom through the leading perturbation if their angular momentum is smaller 
than the one of the perturbation. On the other hand, the KK modes with greater angular 
momentum can cascade into KK modes with lower angular momentum already at zeroth 
order in the isometry breaking, if kinematically possible. Wether all the long-lived modes 
can decay this way is clearly very sensitive to the details of the actual KK spectrum. The 
remaining long-lived modes, if any, have to decay through a sub-leading isometry breaking 
perturbation, or at higher order in perturbation theory, resulting in an exponentially longer 
decay time. 

We note also that, at least if one considers an D3-brane, an isometry breaking perturba- 



2l|, and 



tion w of the warp factor acts as a potential for the angular moduli of the brane 
therefore may determine its position in the internal space. The isometrics left invariant by 
the brane then depends on the perturbation considered, and in particular they may corre- 
spond to isometrics left invariant also by the perturbation. Suppose for instance that the 
leading isometry breaking perturbation has principal quantum number Z = 1 on a 2-sphere, 
^ ~ Y1\4=-i^m'^i^ ■ It is invariant under rotations around a certain axis. The Z)3-brane 
may be driven to the minimum of the potential w, which in the present case (i.e. for I = 1) 
is located on the same axis. Then both the brane and the perturbation leave the system 
invariant under rotations around this axis. In such a case, the long-lived modes have again 
to decay through another isometry breaking perturbation. 

If the throat is approximated by a patch of AdS^xT^'^ , then the leading isometry breaking 
perturbation is a tensor mode 6fij on T^'^ with a ^ 1.29, and potentially dangerous spin 
2 KK relics are listed in rows 2 to 4 and 6 to 11 of Table I. Which of these are long-lived 
depends on the way the SM brane degrees of freedom are embedded into T^'^. The modes 
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with m < 2mo (i.e. those with < 2 x 3.83 in Table I) cannot decay directly into two 
other spin 2 KK modes. Their decay into SM brane degrees of freedom through the leading 
isometry breaking perturbation may be allowed by the last term in Eq. (HTj) . Otherwise, 
they have to decay through subleading isometry breaking perturbations. For instance, all 
the modes of Table I can decay into brane degrees of freedom through the perturbations 



f^^^fki (scalar modes on T^'^) listed in Table II of 2l|], which preserve Lorentz invariance 
and supersymmetry, because these perturbations carry the same quantum numbers as the 
KK modes in question. In this case, 

VXnL ^ SO that the perturbed wave-function of 
the modes has a non-trivial component along Qo = constant, which never vanishes at the 
position of the brane. 

The coupling constant fl551) for the decay of a massive KK mode into two lighter brane 
degrees of freedom through an isometry breaking perturbation with radial profile e~"^/^ 
may be calculated with (H9l) (see also the comment below Eq. (l48l) ) 

Ai^Lb - (^) e'y^'^JAmn'L'Rey^'^). (65) 

In principle, we should sum over all the terms (n', L', M') with non- vanishing matrix element 
in flB5]) . but the term with mass mn'L' closer to the mass rrtnL of the decaying KK mode 
dominate. Again, (|65l) reduces to ( l62l) when the brane is at the tip of the throat, yb = Ut- 

To summarize, we argue that the leading decay channels for the long-lived modes are the 
ones into SM brane degrees of freedom. However, the decay of different KK modes may have 
to be induced by different isometry breaking perturbations of the throat. Different isometry 
breaking perturbations, in turn, lead to different values of a in fl65l) . Note also that throat 
models with different isometrics than the KS throat would lead to different values of a. For 
these reasons in the following we will keep a as a free parameter. The value a ^ 1.29 for the 
leading perturbation considered in \21\ and related to scalar metric perturbations (shown as 
the vertical, "background" leg at the right panel of Fig [2]) will serve as a reference value. 
We will denote by Ag*^^ the coupling constant fl65l) when the brane is at the tip of the throat 
{Vh = yt) 

This is the main formula that we will use to estimate astrophysical effects of the KK 
modes decay. Remarkably, the coupling (166!) can be significantly weaker than the gravita- 
tional coupling 1/Mpi, which makes its astrophysical applications much more interesting. 
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Indeed, while the factor {^) is greater than one, the exponential factor e''^""^^*/'^ can 
be exponentially smaller than one. Already for a throat with warp factor 10~^ and the 
modes with a = 1.29, this factor is about 0.1. Increase of a further decreases this factor, 
for example, for a = 2.29, it is about 10~^. On the other hand, a displacement of the brane 
away from the tip increases this factor. 



V. RELIC ABUNDANCES 



In this Section we estimate the relic abundance of the KK modes. Reheatinsr after 
b..a,Wa,>t,-b..a„e inflat.ou proceed, th.oogh seve... stages Q. The b.a„e/a„ti-b.a„e sys- 
tem mostly annihilates into heavy closed strings. The closed strings very quickly decay 
into the lighter KK modes Q, [isl . The KK modes are significantly lighter than the closed 
strings, and hence they are relativistic when they are produced. They then interact among 
themselves and with the (lighter) SM degrees of freedom on the brane. We can neglect the 
small isometry breaking of the throat for such considerations. Then, the only KK modes 
directly coupled to brane fields are those with zero angular momentum along the directions 
whose isometrics are left unbroken by the brane. For brevity, we denote such modes by KKq, 
and the other KK ones (those that are not directly coupled to the brane) by KK^,. Due to 
their different behavior, the two types of modes are studied separately in two Subsections 
ITOl and IVB] below. 

From these interactions, one can obtain the relaxation time needed for the KK modes to 
.each thecal e,un,b.„.. It has bee. est,„,ated in UQ that. fo. a short (,.flatio„a,v) 

throat, the relaxation time is quicker than the decay of the KK modes into SM fields. The 
situation is more model dependent in the long throat case; in the example of a long throat 
considered here, we fix the parameter A/" by requiring that the warp factor "generates" the 
electroweak scale at the tip of the throat, cf. Eq. (IT^ . However, our discussion below 
is not constrained by any particular inflationary mechanism. For instance, it is possible 
that the SM fields are in the long throat, while the inflationary brane/ anti-brane system is 
in another throat. In this case, the KK modes and the SM fields in the long throat are in 
thermal equilibrium only if the tunneling or the decay of the KK modes from the inflationary 
to the long throat is fast enough; this issue is still under investigation (notice, for example, 
that different conclusions on the transfer of this energy density have been reached in 15 1 
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and [28|). We will briefly discuss these issues in the last Section. 

Our goal is to give a general discussion of the phenomenological limits due to the long 
lived KKr modes. To keep the discussion as general as possible, motivated by the flndings 



of 



15j, in the following analysis we simply assume that the KK modes are in thermal 



equilibrium among themselves and with the SM flelds at some given temperature T. In 
this way, their number density is independent on the history preceding the establishment 
of the thermal equilibrium, and it is simply set by the freeze-out temperature at which 
the interactions changing their number go out of equilibrium. We show below that, in the 
region of parameter space of relevance for phenomenology, this freeze out temperature is 
typically signiflcantly lower than the mass of the KK modes, so that their abundance is 
strongly Boltzmann suppressed. As we mentioned, it is an interesting and rather involved 
problem to estimate the initial number density of KK modes (produced by cascading decays 
of excited closed strings), and to figure out how far from thermal equilibrium the KK modes 
will be. If complete thermalization is not achieved, one should expect a higher number 
density for the KK modes, since the Boltzmann suppression will be reduced. For this 
reason, the phenomenological bounds that we will discuss below based on the assumption 
of initial thermal equilibrium should be considered as conservative ones. 



KKq modes 



Let us first study the KKq modes. We consider both three- and four-legs interactions 
with the brane fields. The trilinear interactions have the coupling A3 discussed in Eq. fl58l) : 
it leads to the decay rate to gauge bosons V, (Dirac) fermions ip, and the Higgs scalar H on 



the brane 



KKo^VV 



A3 



These expressions assume that the KK modes are much heavier than the particle they decay 
into, which is always a good approximations for the cases we are considering (see below). 
The complete expressions, with the kinematical factors included, can be found in 421]. The 
total decay rate is 



A3 



KKo^HH 



A3 



(67) 



KKo^6b 



12 



3x6 + 3 + 3x1/2 1 



160 TT 



320 TT 



+ 



960 TT 



0.047 A^m^ , 



(68) 
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where b stands for the SM particles at the brane. For the inverse decay, we instead estimate 



nb^KKo-W-'XiNb, (69) 

where Nh ~ 0.1 T'^ is the number density of any relativistic species on the brane. This 
estimate is actually valid for temperatures T comparable or greater than the mass m of 
the KK modes. At lower temperatures, the rate has an additional exp {~m/T) suppression 
since most of the light degrees of freedom are not energetic enough to create a KKq mode 



(see for instance the analogous process discussed in appendix B of |39l]). Assuming that 
the energy density of the universe is dominated by the brane degrees of freedom (as we will 
see, this is certainly true at temperatures comparable and lower than the masses of the KK 
modes), we have the expansion rate 

1/2 3 0^/2 

H = , (70) 

where g^, is the number of light degrees of freedom on the brane (we take (7* ~ 100 , as it is 
the case for the standard model at high temperatures). 
This leads to 

^-^^[t) ^^^^P ' 

r^fe^KKo _ 2^^^ ^ T>m. (71) 



H V m 

where 



A^mMp~ 220Ar ( ^V-^ ) . (72) 



5 

We see that both the decay and the inverse decay are effective (namely, F > H) when 
T ^ m if R/ ^/a' ^ 11 for the short throat case, and R/ ^Ja' ^ 250 for the long throat case. 

The four-legs interactions are proportional to the coupling constant A4 = A3 given in 
Eq. (!60|) . It is instructive to compare the rates of the four- and three-legs interactions. For 



To obtain this quantity we have first used the expression ([58|) for A3, ([30| for the mass m of KK modes, 
and (|13p for the ratio Vg /M^ . The remaining parameters combine to give one power of A/" , as defined in 
eq. p6| . The coefficient , has been normahzed to the typical value assumed by the lightest KK modes 
in the spectrum. 
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T ^ m, we find 



TkKo KKo^bb 

rRKo^bb 
Tfe b^KKo KKq 
^bb^KKo 



10-^ XjT^NKKo 9* 
10-2 A2m3 

10-2 XjT'N, 



0.5 A^m^ (- 



5 



T > m , 




(-y ' 



T > m , 



(73) 



where 




2 



8 



A^m^ ~ 5 X 10^ 



(74) 



In general, the 2^2 scatterings dominate at high temperatures, and the moment at 
which the trihnear interactions become dominant depends on the value of R/^/a' . One can 
also check that, at the temperature T c:^ m, these rates are also greater than the expansion 
rate flTOj) . for R/y/a' ^ 6 in the short throat case, and for Rjyfa' <, 37 in the long throat 
case. If this is the case, such interactions lead to thermal equilibrium between the KKq 
modes and the brane fields. 

As the temperature drops below the mass of a KKq mode, both the ratios (173|) acquire an 
exp (— m/T) suppression factor (in this first case this is due to the Boltzmann suppression 
of iVKKo) ill the second case this is due to the fact that the interaction at the numerator 
produces one more heavy particle than the one at the denominator). Also the inverse decay 
is exponentially suppressed in this regime, as we have discussed above. Therefore, the only 
interaction that remains active in the nonrelativistic regime is the decay of the KKq, which 
quickly eliminates these modes. 

B. KKi modes 

The decay rate of the angular KK^, modes into brane fields is suppressed by the small 
isometry breaking of the throat. The decay rate is still given by Eq. fl68l) . but now the 
coupling A3 is given by Eq. (jHS]). Assuming that the brane is at the tip of the throat (?/b = yt), 
and using the relations ffT^ . ffTBl) . and fl5U]) . we can rewrite the total rate normalized by the 
KK mass as 



For the 2^2 interactions, we assume that the s-wave process is unsuppressed; for the present scattering, 
we have actuaUy computed the cross sections for the transverse-traceless hehcity of the KKq modes. The 
cross section for the other hehcities is increased by additional powers of (T/m) in the T > m regime. 




(75) 
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Let us now estimate the relic abundance of these modes. Under the assumption of 
isometric throat, the KK^ modes are not coupled to the brane fields. However, they interact 
with the KKq modes studied in the previous subsection. As long as the modes are relativistic, 

we have 

Tkk-vKkkk ~ 10-2 A^m^ 



T >m : { Tkkkk^kk ~ IO-^AI^Vkk (76) 

i KKKK-^KKKK 

We do not need to distinguish between KKq and KK^ modes in these estimates. Anal- 
ogously to what we saw in the previous subsection, the 2^2 scatterings are typically 
the dominant interactions in the relativistic regime. Also in this case, the rate for such 
interactions is faster than the Hubble rate. Indeed 

rKKKK^KKKK 2 X 10^ f ^ V f ^ V^V- V , T > m , (77) 



H V 5 / V aAV / \m 



which is greater than one for T > m and for R/\'a' <, 4.6 in the short throat case, and for 
R/\^ <, 26 in the long throat case. In the following, we assume that this is the case, so 
that also the KK ^ modes are in thermal equilibrium with the brane fields in this regime. 

These reactions slow down and eventually freeze out when the KK^ modes become non- 
relativistic. In this regime there are three relevant interactions that can potentially reduce 
the abundance of a KKl mode. The first is a 2 — > 2 scattering with brane degrees of freedom, 
mediated by a virtual KKq mode. The second is a 2 ^ 2 scattering into two KKq modes. 
The third is a 2 — > 1 inverse decay, producing a KKq mode. Once formed, the KKq modes 
then quickly decay into the brane degrees of freedom (as we saw in the previous subsection). 
In all these cases, the two incoming KKl modes must combine to produce a zero angular 
momentum (in the directions corresponding to the isometrics unbroken by the brane) . This 
is also the reason why a single KK^, mode cannot interact with KKq or the brane modes. 



Also in this case, the estimate for the 2 — > 2 scattering is given for the transverse-transpose helicity; 
scatterings between different helicities have higher powers oiT/m. 
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The rate for these interactions can be roughly estimated as 



T<m: <( rKKiKK_i^KKoKKo ~ 10"^ Al A^KK^m^ (78) 

Trkl kk_l^kko ~ 10^^ A| A^kkl 5 if sufficiently high momentum 

Since A4 ~ A3 , the first interaction dominates over the second one, due to the greater 
amount of available channels (this is certainly true for the lightest KK^ modes, so that there 
are few lighter KKq modes they can annihilate into; as we will see, the abundance of the 
heavier KK/^ modes is anyhow suppressed relative to the lighter ones). The last process can 
in principle be faster but (due to kinematical reasons) it typically involves only a negligible 
amount of KK^ particles. The two incoming KK^ particles can inverse decay only into a 
heavier KKq species. More precisely, denoting by mo the mass of the KKq mode, and by m 
the one of the incoming KK^, particles of a given species, the momentum of the two particles 
in the center of mass frame is forced to be 



fc=y!^-m2. (79) 

As we saw, the masses of the KK modes are discretized in units of 

Am ~ ^e-^*/-^ .. (80) 

It may be possible that, for some special case, there is a KKq mode having a mass mo slightly 
greater but very close to 2 m . In general, however, the value in ( 1791) will be of the order 
of Am. For the lightest among the KK^, modes. Am ~ m. In the nonrelativistic regime, 
the abundance of the KK^ mode that we are considering is Boltzmann suppressed. Most 
of such particles have momentum ~ T ^ m. The inverse decay KK/,KK_i KKq can 
only happen if, in the center of mass frame, the particles have the much higher momentum 
f l79|) . So, only a very few of the KK^ modes will be able to perform the inverse decay. The 
remaining particles do not have enough energy, and can be depleted only through the first 
two interactions in ( ITHI) . 

For this reason, we neglect the effect of the inverse decay in the abundance of the lightest 
KK^ modes (in this, we differ from 



15|, where the suppression that we have just discussed 
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was not considered). We therefore need to estimate the moment at which the first process 
in flTSl) freezes out. From this, we then obtain the rehc abundance of the KK^, modes. In 
the nonrelativistic regime, 



3/2 



2tt 



-m/T 



(81) 



Using the expression (ffOl) for the Hubble rate we find (see the previous subsection for more 
details in obtaining the final estimate) 



10^ AT ( ^ 

5 



15 



'a' 



R 



m\i/2 
TJ 



-m/T 



(82) 



Equalizing this to unity, on can estimate the freeze-out temperature Tf. 

In Figure [3] we show the ratio of m/T when the scatterings go out of equilibrium in the 
case of parameters which are typical of either a short or a long throat. We also set a rough 
bound of m/T ^ 3 for the nonrelativistic approximation for the KK modes to be valid [40 1. 
For higher values of R/ \fa' the KK^ decouple while still relativistic, which results in a 
greater abundance (at even higher values of R/^foi' we could have that the scattering of the 
KK modes are never in thermal equilibrium; the abundance of the KK particles then depends 
on the specific details of the reheating scenario; due to this strong model dependence, we do 
not consider such cases here). 
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FIG. 3: Estimate for the m/T ratio at which the scatterings KKj^ + KK_i ^h + h freeze out. The 
value is strongly sensitive to the ratio R/^foJ . 
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Since this is the fastest interaction decreasing the number of KK^,, the moment at which 
this interaction freezes out sets the rehc abundance of this species. In general, the later 
the freeze out occurs, the lower is the relic abundance of the mode (since the Boltzmann 
suppression is greater). From the horizontal axis of the Figure [3] we see that, for any fixed 
value of Rj \fa' ^ the higher the mass of a KKq mode (greater ^nh)-, the higher is the value 
of m/T at which the freeze out occurs (later time). Therefore, the KK^, modes with higher 
mass end up with a lower relic abundance. We show this in Figure HI where we plot the 
ratio Ikk^ = ^kk^/s between the resulting number density A^kk^ and the entropy of the 
relativistic fields on the brane 
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(83) 



for the two cases t,nL = 5, 10 . Also in this plot, we show only the cases for which the KK^, 
decouple while nonrelativistic. 
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FIG. 4: Relic abundance for the KK^ species as function of parameter R/^/a' . We note that the 
lightest modes are more abundant than the heavier ones. 



We can verify that the energy density of the KK^ modes is negligible with respect to that 
of the brane fields at the time at which they decouple (we used this assumption in setting 
the Hubble rate (!70l) ). Since the KK^, are nonrelativistic, their energy density is simply their 
number density times their mass. For SM radiation (i.e. the light degrees of freedom on the 
brane), the energy density is instead related to the entropy density through p = (3T/4)s. 
This gives 

(84) 
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It is also worth noting that, for the specific case of warped brane/anti-brane inflation, 
the temperature after inflation is never higher than the mass of the KK modes. Indeed, 
even assuming that the thermalization occurs on a very quick timescale, the energy density 
during inflation, p ^ e'^^^^^ / [A-K^Qga''^)^ leads to a maximal reheating temperature T^ax ~ 
0.23e^^'/'^/-\/a'. Using the expression (l30l) for the KK mass, we have 



T 

J. rf 




(85) 



As we can see in the Fig. HI this value is comparable with the freeze-out temperature that 
we have obtained above. Therefore, our estimate for the number density is reliable in this 
case only if the thermalization is quick, as it appears from the estimates of the relaxation 

nn 

time given in [1^, 115| . 

VI. PHENOMENOLOGICAL CONSTRAINTS 

We are interested in the cosmological effects of the KK^ modes. Such modes are long- 
lived, since their decay rate is proportional to the small isometry breaking of the throat. We 
start by expressing the mass of the KK modes m„ in a more explicit form. By combining 
the three relations f|T3|) . f|T6l) . and fl30|) . we can write 

As we have seen in the previous Section, the lightest KK modes are in general those with the 
greatest abundance. Therefore, for definiteness, in the remainder of this analysis we consider 
only the presence of a single species with C,nL = 5 , as it is typical for the lightest modes 
(stronger effects, and more stringent bounds, are obtained if more KK^ are generated with 
a comparable abundance). The typical values for A/" are given in Eqs. ( fTTll and (fTSl) . In the 
range oi R/ \fa' considered in figures [3] and H] above, we find 

m ~ (2 X 10^^ - 10^^) GeV , short throat 

m ~ (200 - 5, 000) GeV , long throat (87) 

As Eq. f lH^ shows, the KK^ particles dominate the energy density of the universe at 
sufficiently low temperature, provided they have not decayed before that. From the value 
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of Ykkl shown in Fig HI and from the expression (130|) for the mass of the KK particles, we 
can find the temperature of the SM radiation, denoted by Tdom , at which the KK particles 
start to dominate^^. The value Tdom as function of the parameters is shown in Figure for 
comparison, we also show the value of the CMB temperature Teq ~ 7.4 • 10^^° GeV when 
matter starts to dominate in the universe 
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FIG. 5: Temperature of the universe at which the reUc KK^, particles start to dominate, (provided 
they have not decayed yet) as function of R/Va'. The temperature T^q ~ 7.4 • 10^^'^ GeV at 
the moment of matter-radiation equality is also shown for comparison. In the long throat, for 
Rl\fa' ~ 6 , the KK modes can be the dark matter candidate. 




Notice that intersection of two curves, T^q and Tdom for the long throat selects the value 
of the parameter Kj \fa' ^ 6, in which case the KK modes can be the dark matter candidate. 
As we wiU see below, this also poses restrictions on the possible values of a . The Figure [5] 
will be discussed in details in the next two Subsections, where the two types of throats are 
studied separately. 



There is actually a slight difference in the value of entering in the formula for the energy density and 
the entropy after annihilation [i^. This introduces a factor 1.16 in Eq. ((84l) for T <, 0.5 MeV. This 
effect is negligible at the level of accuracy of the present estimates. 

Cosmological quantities in our universe are obtained by using the WMAP values — 0.73, — 
0.27, /i~ 0.70 [41|. 
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A. Short throat 



It is clear from Fig. [5] that, for the case of a short throat, it is phenomenologically required 
that all the KK^ modes decay. If this not the case, our universe would have become (and 
remained) matter dominated much before Zeq ~ 10^, the redshift of matter-radiation equality, 
indicated by the observed values of Qm and of the CMB temperature (in other words, 
today would be too big). 

We therefore assume that the KK^ particles decay with the rate (!75|) . We treat a as 
a free parameter. Typically, smaller values of a correspond to a greater break of isometry 
of the throat, and, consequently, to a quicker decay. From Eqs. (175]) and we find the 
lifetime of the angular KK modes 




We see that this lifetime is strongly dependent on the parameters of the model. For 
illustration, we show it in Figure [6] for a fixed value of / R. 
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FIG. 6: Curves of lifetime (labelled by seconds) of the KK^, particles in the short throat case, as 
function of a and Rj^foJ . For definiteness, we have fixed V^^'° jR = 1 in this plot. Higher values 
of this ratio correspond to a shorter lifetime, cf. Eq. (j88|) . 

Unless a is very high (very long lifetimes), the strongest phenomenological limits are 
those related to the effect that the decay products have on the abundance of light elements 
formed at Big-Bang Nucleosynthesis (BBN). Different limits are obtained for radiative and 
hadronic decays. KK^ have both type of decays, with O (1) branching ratios [42,]. If we 
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separate in Eq. (l68l) the rates into quarks and gluons from the other particles, we find a 
hadronic branching ratio of 0.73 . Ahhough the radiative decays are those which have been 
traditionally more studied, hadronic decays lead to stronger bounds. The BBN limits for 
a pure hadronic decay can be found in the Figures 38, 39, and 40 of |29]. These three 
figures present the constraints for decaying particles of mass m = 100, 1,000, 10,000GeV, 
respectively, in the t — rriY plane (for any given lifetime r, there is an upper limit on the 
allowed values mY). The three cases do not show substantial differences from each other 
(at least, not at the level of accuracy of the present computation). 



The limits of |29|] have been approximately reproduced in our Fig. [10] below. We notice 
that the BBN limit is strongly relaxed and eventually disappears for lifetimes shorter than 
~ 10~^ s : if a particles decays before this time, the decay products thermalize before BBN 
starts. In our case, the product mY is an increasing function of R/^/a', while it is inde- 
pendent of the other "free parameters" of the model {Vq^^/R and a). For the values of 
Rl\fa' we are considering, mY ranges from ~ 400 GeV to ~ 3 x 10^^ GeV. Such values 
are sufficiently high so that consistency with BBN requires that the decay took place place 
before ~ 10"^ s 



It should be noted that the results of [29j have been obtained for much smaller masses than the typical 
KK masses for the short throat; nonetheless, we see that for all the masses considered in [29| (cf. their 
Figures 38, 39, and 40), the BBN bound disappear at r <, 10^^ s . For this reason, we take such value as 
a rough bound in the case of the higher injection energy considered here. 
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FIG. 7: BBN limit r < 10~^ s for three reference values of the ratio Vq^^ /R for the short throat. 
For each case, regions above the lines (greater lifetime) are excluded. The horizontal line at 
a = 1.29 corresponds to the reference value. 

In Figure [7] we show the values of parameters leading to r = 10^^ s . The three lines refer 

1 /6 

two three reference values of the ratio Vq / R. For each case, values above the lines shown 
result in a greater lifetime, and are therefore excluded. 

It is possible that the KK^, particles dominate before decaying. This gives rise to an 
intermediate stage of matter domination. As long as this stage takes place before BBN, 
such a possibility is phenomenologically allowed To see whether this is the case, we 
compute the temperature of the thermal bath on the brane when the decay takes place; by 
equating the decay rate fl75l) (with the expression fl86|l for the mass) to the Hubble rate fITOl) . 
we find 

r,„., ^ 2.2 X lO- GeV (4 x 10') (^) (^^^ (^f^ " . (89) 

In Fig. [5] above we showed the temperature T^om of the brane at which the KK^, particles 
dominate, if they have not decayed yet. Therefore if the value of Tdecay that we have just 
found is greater than the one of Tdom shown in that figure, the KK/, modes decay before 
dominating; the decay increases the temperature of the thermal bath by a negligible amount. 

Such stage can have interesting consequences for cosmology; for instance, it affects the relation between 
the wavelengths of observable cosmological perturbations and the moment when this mode leaves the 
horizon (cfolds number) during inflation. 
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If, on the contrary, Tdom > T^ccay, the KKx, dominate the energy density of the universe when 
they decay. In this case, Tdecay actually refers to the temperature of the thermal bath formed 
by the decay products which dominates over the pre-existing one on the brane. 
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FIG. 8: Different scenarios for the decay of the KK^ particles in the short throat, for the fixed 

1 /6 

value Vq /R = 1 . The horizontal line corresponds to the reference value a = 1.29. 

We show this effect in Figure [HI for the specific case of Vq^^ / R = 1 . The upper line is 
also shown in Fig. [71 Points above this line lead to a lifetime > 10~^ s , and conflict with 
BBN bounds. For points in the intermediate region, the KK^, particles dominate the energy 
density of the universe when the decay. Prior to the decay, the universe is therefore matter 
dominated. 

B. Long throat 

Proceeding as in the previous Subsection, we now find the lifetime for the KK^ particles 

r ^ ^ 3 X 10" s (8 x w^^Y^^'^-''^ ( V"'" ( R 

The lifetime is even more sensitive to the parameters of the model (particularly, to a) 
than in the previous case. For illustration, we show it in Figure [9] for the specific value 

1 /6 

Vg /-R = 10 (as clear from Eq. (1901) . the lifetime is a decreasing function of this quantity). 

The thermalization of the decay products takes place on a much quicker timescale than the decay of the 
KKl particles, and can therefore be considered as instantaneous for this discussion. 



BBN excluded 




Decay while dominant 
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We see that the hfetime changes by several order of magnitudes in the range of value for a 
shown (notice that this range is smaller than the one shown in Fig. [6] for the short throat 
case) . 
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FIG. 9: Lifetime (in seconds) of the KK^, particles in the long throat. For definiteness, we have 
fixed Y^^'° jR = 10 in this plot. Higher (shorter) values of this ratio correspond to a shorter (longer) 
lifetime, cf. Eq. pi|) . 

We also observe that the lifetime of the modes can be easily comparable or much greater 
than the present age of the universe (to — 4.3 x 10^"^ s). The value of a at which the 
two timescales are equal is only weakly (logarithmically) dependent on the other two "free 

1 /6 

parameters" of the model; for instance we find a ~ 1.4 for ^6 = 1. and for aU the 
values of R/ y/a' considered in the long throat case; this is not much different than the value 

1 /6 

a ~ 1.8 that can be seen in Fig. [9] for the Vq / R = 10 case. Due to the strong dependence 
of the lifetime on a, slightly higher values of this parameter result in a lifetime which is 
much longer than the present age of the universe. 

In Figure [10] we show the phenomenological limits on an unstable particle of mass in the 
range (187|) for the long throat case. Limits from BBN are taken from [291 (see the discussion 
in the previous subsection), while the ones from the diffuse 7 ray background are taken from 
Fig. 12 of 3l|- In both cases, the strongest limit comes from hadronic decays, which, for 
the KK^, particles, have a branching ratio of ~ 0.73 . In both cases, the phenomenological 
limits show a weak dependence on the mass of the particle in the interval of our interest; 
this dependence is negligible at the level of accuracy of the present analysis. We also show 
the limit obtained by requiring that the energy density of the KK l particles is smaller than 
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that of dark matter in our universe (this hmit is relevant only for lifetimes greater than the 
present age of the universe). 

Such limits have to be compared with the values of r and of mY obtained for our model. 
For illustrative purposes, we also show in Fig. [TU] the values obtained in the specific case of 

1 /6 

^7i? = 10 (same choice as in Fig. Each line shown corresponds to a given value of a, 
while R/y/Ty' ranges in the interval [1, 29] (that is, the usual one for the long throat case). 
As we already mentioned, the lifetime strongly increases with a. 



1e-05 



> 

S Ie-10 
E 



1e-15 



1e-20 



BBN 


y ray 
















1 .29 ; 


1.6 


i 1.9 



100000 1e+10 1e+15 1e+20 1e+25 1e+30 
T [sec] 

FIG. 10: Combined upper bound on the mass m times the abundance Y for an unstable particle 
with lifetime The BBN limit is taken from 29(], while the bound from the diffuse 7 ray back- 



ground from [31]. The horizontal line denoted by r^mo is the limit imposed by requiring that the 
energy density of the KK particle does not exceed the dark matter 17. The other curves represent 
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the values of r and mY obtained in the long throat, for jK = 10, for different values of a 
(indicated on the lines) and of R/yfa' {jnY is an increasing function of R/Vc/). 

The quantity m Y only depends on R/ yfa' , and it is a growing function of this ratio. 

1 /6 

Therefore, for any fixed values of the other "free parameters" of the model (VJ3 jR and 
a), only values of R/\fa' smaller than some given quantity are phenomenologically accept- 
able. We show this in Fig. [TTl The three lines shown correspond to three reference values 

1 /6 

of V^'^/R. For ah the cases, the points on the right of the corresponding line are phe- 
nomenologically excluded. For values of R/\fa' saturating this bound, the KKj;^ particles 
constitute the dark matter of our universe. In each case, the limit on R/\fa' is strongest 
for an intermediate value of a; this is the value for which the lifetime of the KK l particles 
is comparable to the present age of the universe (we see in Fig. [10] that this is the point at 
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which is most constrained). At higher values of a, the hfetime of the KKl particles 
becomes much greater than the present age of the universe. In this case the only relevant 
bound is that their energy density does not exceed the one of dark matter in the universe. 
Since the energy density depends only on R/\^, we find the same limit Rj \foi' ^ 6 for all 
the values of a and V^^^ / R in this range. When this limit is saturated, the KK modes are 
identified with the dark matter of the universe. On the contrary, for the lowest values of a 
shown, the modes decay before the onset of BBN, and the phenomenological limit disappear. 

1 /6 I 

The "scaling" of the curves with Vq / R shown in Figure [TT] can be also easily understood. 

_____ /q 

From Eq. fl^Ul) . we see that the lifetime is a decreasing functions of / R. This decrease 
is "compensated" by the growth of a in the numerical prefactor (which is by far the most 

1 /6 

sensitive quantity on a in that relation). Therefore, as VJj /R increases, the same lifetime 
and the same "pattern" in the phenomenological limit is obtained by a small (logarithmic) 
increase of a). We see, in general, that for any given value of a, the decay can take place 
before BBN, provided the volume of the compact space is sufficiently large. 
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FIG. 11: Final exclusion region for the parameter in the long throat. The three lines correspond 

1 /6 

to three reference value of Vq /R. For each case, values of the parameters on the right of the 
corresponding curve conflict with the phenomenological limits shown in Fig. [TOj The highest values 
of a shown result in KK^ particles with a much longer lifetime than the age of the universe. In 
this case, the only relevant bound is that the energy density of the KK^ particles does not exceed 
the one of dark matter in our universe. The horizontal line corresponds to a = 1.29 

We conclude by commenting on the case in which the angular KK modes are identified 
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FIG. 12: Angular KK modes can be the dark matter candidate for Rj\foL' ~ 6 and for the values 
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of the parameters a and Vj /R above the line shown in the plot, corresponding to a lifetime of 
the KK particles greater than about 3- 10^^ s (to avoid the limit from the diffuse 7 ray background, 
cf. Fig.Hni). For V^'^/R > 1 this occurs for relatively large values of q > 1.7. 

with the dark matter of the universe. As we have aheady mentioned, the right abundance 
is obtained for R/ y/a' ^ 6 . Moreover, hmit from the diffuse 7 ray background require that 
the hfetime of the modes is greater than about 3 ■ 10^^ s (cf. Fig. [TOl) . The hne in Figure fT2] 

1 /6 

represents the point in the Vq /R — a plane corresponding to this lifetime. Points below 
this line give a shorter lifetime, and are therefore phenomenologically excluded. In particular 

1 /6 

we see that a > 1.7 is required for any value of Vg /R > 1. Therefore, the reference value 
a = 1.29 does not allow to identify the angular KK modes with the dark matter of our 
universe. 

VII. DISCUSSION AND CONCLUSION 

In this paper, we have studied the problem of angular KK relics in string theory cosmol- 
ogy. As a specific example we considered angular KK modes resulting from approximate 
isometrics of the internal space in the context of warped flux compactifications. We focused 
on the case where the isometrics of the warped throat are distorted by its embedding into the 
compact manifold. The profile of the different isometry breaking perturbations of the throat 
(the parameters a/, in (fT2l) ) has been studied in ^] for a Klebanov-Strassler background, 
and we took their results as reference values. 

We investigated the interactions and possible decay channels of the KK modes with an- 
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gular momentum, calculating the corresponding coupling constants. We did explicit calcu- 
lations for the spin 2 KK modes. We found that the lightest modes with angular momentum 
cannot decay directly into any other two KK modes (massive or massless), independently of 
the isometry breaking perturbations. We considered the model where the SM particles are 
localized on a (probe) 3-brane inside the throat. In this model the leading decay channel of 
the angular KK relics is the decay into brane degrees of freedom. Different KK modes may 
decay through isometry breaking perturbations with different exponents e""^*/"^, resulting in 
exponentially different lifetimes. Our final expression for the lifetime of the relic KK modes 
then depends on several parameters: the degree of isometry breaking at the tip of the throat 
^-ayt/R^ the volume of the internal space Vg, the AdS curvature radius R, and the throat 
warping = e"^'/^. For completeness we also included as a parameter the position of the 
SM brane inside the throat e'^^^^. Remarkably, the coupling between angular KK modes 
and the SM particles can be much weaker than the gravitational coupling. 

We then studied the cosmological constraints on the parameters of the single throat 
model. We distinguished two limiting cases: a short throat where brane inflation occurs, and 
a long throat related to the hierarchy problem. In both cases, we considered the Standard 
Model fields to be located at the tip of the throat and to be at some moment in thermal 
equilibrium with the lightest KK modes. The freeze-out abundance of the KK relics then 
depends essentially on a single parameter {R/ v^), although the dependance is very strong. 

We considered the influence of the KK relics on the thermal history of the universe and put 
together observational constrains on the parameters of KK modes and consequently, on the 
underlying throat model, using limits from BBN 29|] and the astrophysical 7-ray background 



3l| , in the case of decaying modes, and limits on the dark matter energy density, in the case 



of lifetimes much greater than the present age of the universe. Such limits are very sensitive 
to the parameters of the models (see Fig. [7] for the short throat, and Fig. [11] for the long 
throat cases), but we can exhibit some general trends. 

In the short throat, the lifetime of the KK relics is much shorter than in the long throat 
(since the warp factor is smaller and, consequently, physical times are quicker), but their 
relic abundance is also much higher. In this case, we must require that these modes decay 
before the onset of BBN, namely with a lifetime r ^ 10~^ s . For greater lifetimes, the 
BBN and overclosure limits are typically exceeded by many orders of magnitude ; therefore, 
the demand that the decay takes place before the onset of BBN holds even if the thermal 
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history for the KK modes is shghtly different than the one considered here (for instance, 
even if they do not reach perfect thermal equihbrium before they decay). We also showed 
that, for natural values of the parameters, the KK modes can dominate the energy density 
of the universe before decaying, see Fig. [HI This leads to an intermediate stage of matter 
domination in the early universe, which may have interesting consequences. For instance, it 
alters the thermal history of the universe, it affects the relation between the wavelength of 
present-day cosmological fluctuations and the corresponding number of e-folds of inflation, 
etc. 

For the long throat case, a large fraction of the parameters space is excluded, see Fig. [TTl 
In this case, for typical values of a obtained in [21] (such as a = 1.29 for their leading 
isometry breaking perturbation), the constraints from the late decays require a sufficiently 
small value of R/^fa' ( 5, or smaller), so that the abundance of the modes is strongly 
(Boltzmann) suppressed (recall that R/ ^/a' ^ 1 is required for the validity of the super- 
gravity treatment of the throat geometry). Otherwise, lower values of a, or a sufficiently 
large compact space {yl^^ jR ^ 100), is required for these modes to decay before BBN. In 
the opposite regime, the lifetime of the KK relics can be much greater than the present age 
of the universe. If this is the case, or if, due to selections rules, some modes are stable, one 
gets the right abundance for these modes to be viable cold dark matter candidates if one 
tunes the value of R/ \fa' ^ 6.1. This also requires a value of a 1.7. 

In our analysis, we considered the Standard Model degrees of freedom to be located on 
a probe 3-brane at the tip of the throat where the KK modes are localized. The ultimate 
fate of the KK relics in a given model will depend on the specific way that the Standard 
Model is realized and embedded into the internal space. Furthermore, different or additional 
constraints may apply if the Standard Model degrees of freedom are located in the unwarped 
region of the Calabi-Yau or in another throat. 

Already a displacement of the SM brane from the tip of the throat would affect the decay 
rate of the KK modes. In Sections 5 and 6, we only considered the probe brane to be at the 
tip, = yt- Otherwise (?/{, < y^, the decay time of the rehc KK modes into SM particles 
would be much bigger, see Eq. fl65l) . Shifting the SM brane from the tip has qualitatively 
similar effects than an increase of a, but quantitatively it is more involved. If y^ < yt, all 
the KK modes interact much more weakly with the brane. In this case, even the modes 
with zero angular momentum can be dangerous relics, and observational constraints would 
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dl+ml^-^)^m = 0, (91) 



be more stringent. 

We restricted ourselves to the single throat case. The reason is related to at least two 
additional complications that arise if we include another throat. One complication arises for 
two throats with very different warpings. The physics of inflation and subsequent thermal- 
ization will corresponds to values of the Hubble parameter and of the temperature which 
exceeds the mass of the KK modes in the long throat (~ TeV). This generates a mass gap 
for the KK modes. The tunneling of KK modes from the short to long throat should take 
this effect into consideration. As it was recently shown in 46||, in such a scenario the Stan- 
dard Model throat (together with SM brane) will be cloaked by a Schwarzschild horizon, 
and the usual treatment of tunneling is not applicable. 

The second complication is even more restrictive. Let us try to model the wave equation 
of the the angular KK modes in the bulk space between throats [l^ . Since warping there is 
insignificant, instead of equation (l23l) we have 

where we continue to use the "radial" coordinate y (in a loose sense) and where is 
the eigenvalue of the Laplace operator Aj in the directions orthogonal to y. Assuming 
that the angular momentum is connected to p^, we conclude that the angular KK modes 
are exponentially suppressed in the bulk, ~ e"^'-^^'^'/'^. For V^^^ / R ^ 1 this effect 
exponentially suppresses the tunneling rate. 

An interesting extension of the present analysis would be a more detailed study of the 
thermal history of the KK modes (or glueballs, in the dual theory), from their production 
by the decay of heavy closed strings formed after brane annihilation, to their thermalization 
among themselves and with light SM fields. In this work, we have conservatively assumed 
that complete thermalization is quickly achieved; this reduces significantly the abundance 
of the KK modes (due to Boltzmann suppression). One should expect the KK modes to 
have a much higher abundance if thermalization is not complete; this would lead to even 
stronger bounds on the parameters of the model than those obtained here. 

We note the paper 47|] , which also made connection between cosmology and the glueballs 
in the dual picture of the gauge theory with extra isometrics, which appeared as the present 
paper was being completed. 
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In this Appendix, we derive the effective action to second and third order in the spin 2 
KK modes ( IT9|) around the general background ([2]), assuming only 4-dimensional Lorentz 
invariance but independently of any isometrics in the internal space. We then show that a 
massive mode cannot decay into one or two zero mode, generalizing a similar discussion in 
Section lYl 

We consider the bulk action 



where the lagrangian L corresponds to all the matter in the bulk (scalar fields, forms, ...), 
as well as possible brane localized sources. If the spacetime has boundaries, we should also 
add the appropriate Gibbons-Hawking boundary terms. The ansatz for the D-dimensional 
metric is 



where capital latin indices M, A, ... run over the D-dimensional coordinates, greek indices 
/i, z/. A, ... run over the 4-dimensional coordinates, and latin indices a, 6, c, ... run over the 
internal [D — 4)-dimensional coordinates of the compact space. 

The background solution with 4-dimensional Lorentz symmetry corresponds to g^^y = rj^y 
in (|93|) . For the spin 2 perturbations, g^j, will also depend on the coordinates and y'^, 



APPENDIX 




(92) 



ds^ = GMNdx^dx^ = e^^^y^^ gt,u{x^, y") dxf'dx" + gab{y^) dy^dy^ , 



(93) 
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while the metric of the internal space gab depends only on the extra-dimensions coordinates 
Einstein equations are 

T SH 

Rmn = Tmn — £) _2 ^^^^ with Tmn = ^ Gain — 2 ^qmn ' ^^^^ 

where T is the trace of the energy- momentum tensor, T = G^^Tmn- The background 
symmetries correspond to a matter lagrangian which is independent of the 4-dimensional 
metric: = 0. The (/i, v) and (a, h) components of Eq. (|9^ then read 

= + e"<,,„, (95) 

and 

«.,._l_(_,+^..^)^„,_2-p,. (96) 

It is straightforward to compute the Ricci tensor of the metric f l93|l . The (/x, z/) compo- 
nents are 

RAG] = <^M-^V,V^(e^V) 

+ \ e-'^ g'^ [2 {e'^g,,) 3^ {e'^g.,) - 3^ [e'^g,,) 3^ [e'^g,.)] (97) 

where R^^} [g] is the 4-dimensional Ricci tensor of the metric g^v, Va is the covariant derivative 
associated with the metric gab-, and a,b,c, ... indices are raised and lowered with the metric 
cjab- Similarly, the (a, b) components of the D-dimensional Ricci tensor read 

Rab[G] = Rab[g] - ^ V. [e-^V^S, (e^V))] - le^'"" g^"" g'' 3, (e^V) db (e^Vp) (98) 

where Rab[g] is the {D — 4)-dimensional Ricci tensor of the metric cjab- To compute the 
effective action, we will also need the D-dimensional Ricci scalar for the ansatz (p3l) . From 
Eqs. fl97l [98]) . and performing the derivatives of the warp factor, we get 



R[G] = [g] + R[g] + i g^'' g^^ {3,g,p 3^g,,^ - 3,g,, 3^g,,) - \ 3^g^^ 3,g,, - 

g^" ^cV'g,u - 3' A {6g^''3,g,, + g,,3,gn - 8 V^V'A - 20 3,A 3' A (99) 

where R^'^^[g] and R[g] are respectively the 4-dimensional and {D — 4)-dimensional Ricci 
scalars of the metrics g^i, and gab- 
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We first consider the background solution, with g^^ = rj^j^^. Eq. (195|) with (l97Il then gives 



D-2 



^-9 



cd 



6C 



cd 



Similarly, contracting Eq. ( 196|) with g and using Eq. ( 198|) gives 



+ 2 £ = 6 VeV'A + 12 



(100) 



:ioi) 



where we have used Eq. fllOOp to eliminate the term in g"''j^. 

We now consider the spin 2 perturbations h^i,{x,y), defined in f|T9l) . To obtain their 
effective action, we start from Eq. (!92l) with (!99l) and v^— G = ^^"^^ may then 

eliminate the background {D — 4)-dimensional Ricci tensor R[g] and matter lagrangian C 
by using the relation (110 ip . This gives 



S — 5*1 + 5*2 



(102) 



with 
and 

^2 ^ 



S, = ^^^ I d^x^g 



-9 e 



2A 



(103) 
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g''" ^c^'g,. - d'A (Gg'^'d^g^, + g^.d^g^ - 2 VS/'A - 8 d,A d'^A 



(104) 



We first expand the action to quadratic order in h^^,. After some integrations by parts 
(and neglecting the boundary terms), we can cast the action in the form 



^ i d'^x^ge^^h^'' 



h 



(105) 



where greek indices are now raised and lowered with the 4-dimensional Minkowski metric 
rj^u and 0(4) = rj^'^d^dv. The vanishing of the terms inside the brackets gives the equations 
of motion to first order in h^y. With the separation of variables fl20l) . this gives Eq. fl22|) for 
the profile of the modes in the internal space. The action fll05p then becomes 



M 



D-2 
D 



8 



E / d^'-Wge'^^m^^' / d'x^^^^-^'^ [U 



J(4) 



m 



(m) 



(106) 



Therefore, with the orthonormal conditions 



d^-'y^ge'^^^^ 



2Ai 



J mm! 1 



(107) 
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the quadratic action decouples, and each 7^™'' represents from the 4- dimensional point of 
view a canonically normalized spin 2 field of mass m. For self-adjoint boundary conditions, 
the orthogonality conditions in fll07p follows from the mode equation fl22|) in the usual way, 
by multiplying ( |22i) by integrating by part over the internal space, and subtracting the 
equation obtained by interchanging m and m' . 

Eq. (I22p admits in particular a massless mode with a constant wave function in the inter- 
nal space, $0 = constant. This mode represents the standard graviton in the 4-dimensional 
picture. The effective 4-dimensional Planck mass may then be obtained by performing the 
y-integral in (1103^ for a 4-dimensional metric Qnuix) depending only on the external coordi- 
nates 

M^i = ME~' J d^SV^ge'^ = ME-' Vn-4 , (108) 
where Vd_4 is the volume of the internal space. Using the normalisation fll07p . the wave 
function of the zero mode is then 

We now expand the action to third order in /i^j,, in order to obtain the trilinear interac- 
tions between the spin 2 KK modes. For concreteness, we consider the decay of a mode of 
mass m into two modes of mass m' and m". The first term in (11021 11031) gives contributions 
of the form 

D ^ I d^-'y^ge'^^A^^,, J d'x^f> dA:r;^d'^ ^^^"^^^ .. (110) 

In general, the y-integrals do not require the conservation of the total KK mass, m = m'+m" , 
as would have been the case in a flat background. However, if one of the decay product is the 
zero-mode, m" = 0, its wave-function is constant and goes out of the integral. The y-integrals 
then reduce to the orthonormal conditions fll07p . which imposes m = m'. Therefore, a 
massive KK mode (m 7^ 0) cannot decay into 2 zero modes gravitons {m = m' = m"), as 
also noticed in ij]. Here we note furthermore that a massive mode cannot decay either 
into another massive mode plus a graviton since, for the process m —>■ m' + 0, the constraint 
m = m' is not compatible with the conservation of the 4-dimensional energy-momentum 
flHT]) . We thus conclude that a massive KK mode cannot decay directly by 3- legs interaction 
into one or two zero-mode(s). 

We can check that the second term in fll02p does not change these conclusions. Expanding 
(I104p to third order in h^i, and performing integrations by part, we find after a tedious 
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calculation 



^^^^ Id^'x^ge'^ 



1 
3 



111) 



We are interested in the contribution of this expression to the process m — m' + 0. There is 
only one possibility for the constant zero-mode in the first term, 2 in the second term and 3 
in the third term. The 4-dimensional coupling constant for this process is then proportional 
to 

i d"^^, - dPA d, $„0 
WcV^A + Ad^Ad^A^ • (112) 

By integrating the first and second terms by part, and by using the mode equation fl22l) . 
this reduces to 



M 



D-2 
D 



(113) 



which again is non- vanishing only for m = m', from the orthogonal conditions ( 1107p . 

Therefore, the 3-legs decay of a massive spin 2 KK mode into 1 or 2 graviton zero-mode(s) 
is not possible. We will use this in Section [IV] when we study the possible decay channels 
of the KK modes. 
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